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ABSTRACT 


This thesis proposes an analytical model to test various 
assumptions about conventional/chemical warfare. emia ' Ss 
status in conventional/chemical combat is modeled as states 
in a semi-Markov chain with transient and absorbing states. 
The effects of differing chemical threat levels, availability 
of decontamination assets and assumed personnel degradation 
rates on expected unit life and capabilities are tested. The 
model results indicate a possible optimal mix of conventional 
and chemical weapons. Also the availability of 
decontamination assets affects expected unit life more than 


decisions as to when to decontaminate a unit. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


Purported Iraqi use of chemical agents on the city of 
Halabja on 4 April 1988 resulted in as many as 5,000 deaths. 
The agents used were reported to be hydrogen cyanide, mustard 
gas and possibly Sarin, a deadly nerve gas [Ref. dE This 
incident highlights the continuing presence of chemical 
weapons in the arsenals of many of the world's armies. Other 
than the Iran-Iraq war and possible use in the Afghanistan 
war, modern armies have not seen chemical warfare on a large 
scale since the end of WW I. With the development of new 
chemical agents in WW II and new weapons to deliver those 
agents the historical data from WW I do not reflect the 
probable outcome of chemical warfare today. 

Information from the Iran-Iraq war serves to dramatize 
the possible effects of large scale chemical warfare. 
However, these attacks against unprotected civilians and 
troops do not directly translate to chemical warfare between 
two modern, fully equipped armies such as the United States 


and Warsaw Pact forces. Limited information from small 


: Hydrogen cyanide(HC) is the gas used in U.S. gas 


chambers to kill convicted criminals. Mustard gas is a 
vesicant which kills tissues on contact, cauSing severe 
blistering of exposed skin and destruction of mucous 
membranes. Sarin is an organophosphate compound which kills 
by disrupting the transmission of nerve impulses. Both HC and 
Sarin can be fatal in seconds to minutes. 


a 


conflicts or scattered use of chemical agents do not enable 
planners to predict the effects of large scale use of chemical 
agents or the tactics that would be needed to minimize the 
effects of such use. 

In planning for the eventuality of chemical warfare the 
tactician's primary tool is simulation and modeling. To be 
useful, these models need to reflect the full extent of our 
chemical warfare knowledge. Such models also need to be 
integrated into current combat models instead of being add-on 
external models. Integration reflects the U.S. Army view that 
any future conflict will be one of conventional and 
nonconventional (nuclear, chemical, biological) weapons. 

The purpose of this thesis is to demonstrate the use of 
an analytical model as an important tool for gauging the 
effects of chemical warfare. The particular model considered 
is a semi-Markov chain with transient and absorbing states. 
The model will be used to determine the effect of the 
following parameters on the survivability and capability of 
a unit in combat: frequency of enemy use of chemical weapons, 
availability/priority of decontamination assets and assumed 


degradation effects on personnel in a chemical environment. 


II. STATEMENT OF THE PROBLEM 


In some current high resolution and aggregated Army combat 
models the effects of chemical agents on the battlefield are 
portrayed, but not necessarily accurately. The parameters of 
interest, such as how the degradation of aeeunit's 
effectiveness due to being placed into Mission Oriented 
Protective Posture (MOPP) affects the outcome and pace of 
battle, are only now being addressed. For most of these 
models, testing the sensitivity of the model to the assumed 
chemical parameters is, at best, very difficult and time 
consuming. Production runs from one of these aggregated 
combat models can take on the order of days to weeks to run 
and the postprocessing necessary to determine the effect of 
model changes on the Measures of Effectiveness (MOES) can take 
weeks to months. Less cumbersome methods are needed to test 
the sensitivity of battle outcome and the pace of battle to 
the use and effects of chemical agents. 

Modeling a unit in combat in a conventional/chemical 
environment as a semi-Markov chain is one way to accomplish 
the needed sensitivity analysis. Using widely accepted high 
resolution input data for conventional conflict, the user can 
then input assumed chemical warfare effects and find a range 


of values that the model is extremely sensitive to or 


extremely robust against. This thesis proposes such a model 
and demonstrates the sensitivity and MOE evaluation process 
using well known techniques involving Markov chains and their 
properties. 

As an analytical model, this semi-Markov model has several 
advantages. Most importantly there is no need for 
replications to reduce variance in the results. Once the 
transition probabilities, transition times and starting state 
are specified, the results of the model (expected unit 
lifetime, probability of absorption, etc.) are 
uniquely determined. Also, the solution techniques for this 
class of model are widely known and the computations are 
relatively fast and simple. However, the output from this 
model must be interpreted with the following caveats. First, 
the model is very scenario specific and may require simulation 
model input for the transition probabilities and times. 
Second, the model is not intended as a replacement for time 
step simulations but rather as a tool for sensitivity analysis 
for parameters that may be later used in a simulation model. 

This thesis will use a semi-Markov model to investigate 
how the frequency of enemy chemical usage, 
availability/priority of decontamination assets and assumed 
personnel degradation rates effect the survivability and 


capability of a unit in a conventional/chemical battle. 


II]. METHODOLOGY 


The following definition of a Markov process is extracted 


from Introduction to Probability Models by Sheldon M. Ross. 


merconstaer a stochastic process ({X_, m= 0,1,2,.22) that 
takes on a finite or countable number of possible values. 
Unless otherwise mentioned, this set of possible values 
of the process will be denoted by the set of nonnegative 
integers (0,1,2,...}. If X, = 7, then the process is said 
to be in state 7 at time nm. We suppose that whenever the 
process is in state 7, there is a fixed probability Ds; 
that it will next be in state j. That is we suppose that 


P(X =JIX= 1, Xpy=Tpeye ese X=T,, Xp=MQ} = Py; (1-1) 
iomaliaestates Wie. ,i api; y sand all n > 0. Such a 
stochastic process is known as a Markov chain. Equation 


(1.1) may be interpreted as stating that, for a Markov 
chain, the conditional distribution of any future state 
feneGaven the past states X,, X,,...,X%,., and the present 
state X,, is independent of the past states and depends 
only on the present state. [ Re figeu2 Digell 3 :2)) 

Note that for a Markov chain the transition times are all 
one time step. The definition above specifically refers to 
a Markov chain with stationary probabilities (i.e., the 
probability of transitioning from state 7 to state j is not a 
function of time). 

For a process described by a Markov chain with stationary 
probabilities, we can answer three questions. First, if there 


are absorbing states (states which the process can transition 


into but can never leave) what is the probability that the 


process is ultimately absorbed in state j given the process 
started in state 7? Second, what is the expected time until 
the process first enters an absorbing state? Third, what is 
the expected number of visits the process makes to each state 
prior to absorption? 

The following method for calculating the quantities of 


interest comes from An Introduction to Stochastic Modeling 
by Taylor and Karlin. This method of calculating the 
absorption probabilities uses matrix algebra and gives an 
intermediate matrix W called the fundamental matrix that can 
be used to determine the mean number of times a state is 
visited before absorption occurs and the mean time spent in 
each state prior to absorption. Die referred to as the nth 
step probability, denotes the probability of being in state 
j after n transitions given that the process started in state 
ee 


Consider a Markov chain whose states are labeled 0, 


lt, eae States 0,1,...,r-1 are transient in that 
Bee) > 0 as n-o for 0 < 7, j < r, while states r,...,N 
are absorbing, or trap, sand hence) etoile) <N eee 
transition matrix has the form 

Q R 

P = (1) 

0 I 
where 0 is an (N - r+ 1) matrix all of whose components 
are zero, I is an (N - r+ 1) X (N - r+ 1) identity 
matrix and q;; = Pi; £ One O<ae icilois [PS <ap ens i o Rema Os Er cy? 


The R matrix in (1) has a row for every transient state 
and a column for every absorbing state in the chain. The 
individual entries, r,;;, are the probabilities of transitioning 
from the transient state 7 into the jth absorbing state in one 
time step and are given by ae eo oe Ost) 2 ee Delng =the 
state label of the first absorbing state as defined in the 
previous quote from Taylor). 

Once the transition matrix is expressed in this format 


the fundamental matrix is defined by: 


W= (I - Q)” (2) 


where I is an (r x r) identity matrix the same size as Q. 

Taylor shows that the values wW,; are the expected number 
of visits to state j before absorption given that the Markov 
chain initially started in state 7 (Ref. 3:pp. 117-118]. If 
the starting state of the chain was a distribution of states, 
then the expected number of visits for any state Jj can be 
calculated by: 


V = WP tial (3) 


where P. 


ee] is a column vector of the probabilities that the 


initial state was state 7? and V is a column vector of the 


number of visits for each state given the initial distribution 
Fiscal 
Now using the matrix W the absorption probabilities are 


calculated by: 


U=WR (4) 


using the R matrix defined in the P matrix in (1). 

The resulting U matrix has a row for every nonabsorbing 
state in the chain and a column for every absorbing state. 
The entries u,, are the probabilities of absorption ingctarme 
J given that the chain started in state 7. [Ref. 3:p. 119] 

All of the above calculations are derived from the pure 
Markov chain but can be used for a more general class of 
stochastic models called a semi-Markov’- chain. The 
distinguishing feature of a semi-Markov chain is that the 
process can sojourn in the state 7 for a random time with a 
mean of pt. before transitioning to some state j [Ref. 2:p. 
292]. Just as with the transition probabilities, the mean 
sojourn time, [,, 1S independent of how the process reached 
the state ?. For the particular model considered in this 
thesis, the sojourn time is a degenerate random variable that 
takes on a constant value with probability one. Transition 
times are assumed to be a function of both the state the 


process is leaving and the state it is transitioning to, being 


constant for given? and j. Thus the uw, for a state 7 is given 


by: 


a P38); (5) 


where the ¢,, is the constant transition time from state 7 to 
state J. The matrix T will denote the transition time matrix 
whose entries are the individual transition times t,;- 

Note that the absorption probabilities (4) and expected 


number of visits per state prior to absorption (1) are 
functions only of the P matrix for the chain. Thus a Markov 
and a semi-Markov process that have the same P matrix, will 
have identical W and U matrices. However, the expected time 
until absorption and the time spent in a state prior to 
absorption will not be equal for the two processes. 

For a Markov process all transition times are one time 
step; therefore, the time spent in any state is exactly equal 
to the number of visits to that state. For an initial state 
7 the expected time until absorption is then the sum across 


the 7th row of the W matrix (6). 


For a semi-Markov chain the expected times until absorption 


are given by: 


Tike =W a (7) 


where #2 1S a column vector of the expected sojourn times for 
each state, W is the fundamental matrix and T,,. is a column 
vector of the expected times to absorption given the process 
starting state. For a semi-Markov chain the expected time 
prior to absorption spent in a particular state j given the 
process started in state 7 is given by W; jh. 

The following numerical example demonstrates the above 
Calculations. We assume a semi-Markov process with five 
states {0,1,2,3,4} where the state space diagram is given by 
Figure 1. The numbered circles represent the states in which 
the process can exist. An arrow from state 7 to state j 
denotes a positive probability of transitioning from state 7 
into state j. The numbers above the arrows are the transition 


probabilities and the numbers below the arrows are the 


transition times. 


Ie 





Figure 1 State Space Diagram for a semi-Markov Process. 


The one-step transition probability matrix for this 


process is given below. 


O t2 2D 220 we 
ORS 18 O 7 S. O 
P = O slo 0295 O 0 
O 0 O 1 0 
0 O O 0 1 


This semi-Markov process also has a transition time matrix 


T shown below. 
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For this example the Q matrix becomes 


0 225 2295 
Q = 5 0 0 
0 ~75 229 


and the R matrix is shown below. 


Using equation (2) the resulting W matrix is shown below. 


W = 0.667 eo 33 PAPA 


1.333 0.667 0.444 
0.667 os <256 


aS 


This "fundamental" matrix is the number of visits prior to 
absorption. The absorption probabilities and the expected 
time to absorption for each initial state are listed in Table 


I below. 


Table I. SAMPLE ABSORPTION TIMES AND PROBABILITIES. 


Probability of absorption 


Tnit ia. ea Expected time 
State State 3 State 4 to absorption 
0 0.667 0.333 Oslo 
al 0.333 0.167 4.389 
2 O3355 OO. 167 8.722 
3 IS 18716) 0 0 
4 6) 1.000 0 


IL 


Note that the U matrix defined by (4) is the first three 


rows of columns 2 and 3 in Table I and the Tp, Vector defined 


Ss 
by equation (7) is the first three entries in the last column 
of Table I. 

The solutions to the example problem and to the full model 


were obtained using the programming language APL on an AST 


personal computer. 


Age 


A. 


BASIC CONCEPT 


IV. MODEL FORMULATION 


In the "real world" a unit in combat is assigned missions 


and succeeds or fails on the basis of the factors given by the 


acronym METiI= ia 


Ibe 


Ze 


3% 


Le 


Dy 


Mission: 


Equipment: 


Troops: 


Terrain: 


Time: 


METT-T is an abbreviation for the following: 


What is the unit's current mission and 
what do we want the unit to do in the 
future?. 


What is the quantity and condition of the 
unit's equipment? Is the on hand 
equipment sufficient for some missions and 
not for others? 


Analogous to the question for equipment, 
what is the quantity and condition of the 
unit's personnel? Is the unit at full 
strength for personnel or at twenty-five 
percent strength for personnel? Included 
in this category are morale, discipline 
and leadership, which are very hard to 
quantify. 


What is the terrain in the unit's area of 
operations and how does this terrain 
affect the unit in the performance of its 
mission? Also if a unit must move to 
another area of operations, how is the new 
or intervening terrain expected to affect 
the unit's operations? Terrain in this 
instance is defined to include the weather 
and its effects on the battle. 


How much time is available for the 
planning and execution of the mission? 
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All of these questions are considered when a real 
commander is either deciding on a mission for his unit or is 
assessing his unit's probability of mission success given a 
specific mission. A commander also considers these factors 
for the enemy he is facing. For a combat scenario, if the 
enemy threat, terrain and overall unit mission remain 
relatively constant, the commander is making decisions based 
on his unit's "state" assuming adequate time. The Ss a Une 
state is determined by the unit's equipment and personnel 
strengths and the unit's current activity. For example, a 
unit that is at 20% strength for both personnel and equipment 
will not be sent back into battle if at all possible but will 
more likely go to the rear for restorative processes such as 
resupply and refit. If, however, the tactical situation 
forces the commander to recommit this unit without 
Bectemationn his unat has a high probability of mission 
icp re . 

Note that in "real world" combat a unit's state expressed 
in terms of personnel and equipment strengths does not present 
the entire picture. A unit that has suffered sustained 
attrition and is at 50% personnel is likely to be more 
effective than a unit that is at 50% personnel due to a 
catastrophic loss. In the former case the unit's chain of 
command has had time to adapt to the loss of personnel while 
in the later case the sudden shock has probably disrupted the 


unit's chain of command and lowered morale more than in the 


25 


former case. However these effects are beyond the scope of 


most models to credibly capture. 


B. SEMI-MARKOV CHAIN APPROACH 

One approach to modeling a unit's transitions from 
state to state is to use a semi-Markov chain. In the "real 
world" units start at various strength levels and cycle 
through processes that lower strength (combat) and processes 
that raise strength (decontamination, resupply/refit). 
Ultimately the battle is over or the units leave the battle 
due to being ineffective or being rotated to another battle 


area. A semi-Markov chain can capture this change in status 


and subsequent absorption. Such an approach assumes the 
following: 
1. The unit can be in specified states and can transition 


between these states. 


2. These unit states can be described in sufficient detail 
so that the next possible states a unit can transition 
into can be determined solely from the unit's current 


state description. This implies that the process is 
memoryless (i.e., how a unit arrived at a state is not 
important). 

3. It is desirable to have these transition probabilities 
remain stationary (i.e., invariant over time) so that 


every time a unit enters a given state its probability 
of entering any other state is the same. 
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Note that the memoryless property fails to consider the 
shock effect, discussed in paragraph A. above. This is a 
shortcoming of the current model but could be overcome by 
expanding the state space to include an indicator that tells 
how a unit arrived at the current state. Units with a 
catastrophic loss indicator would have different transition 


probabilities than units with a gradual loss indicator. 


SaeDEFINITION OF TERMS 

To be useful the state description must be detailed enough 
to allow all of a unit's next transitions to be determined 
based solely on the current state and broad enough so that the 
size of the state space does not become unmanageable. For 
this model, the state description covers three basic areas: 
the current unit activity, the current unit strength and the 
current unit environment. The terms used to describe each of 
these areas are detailed in the following sections. 

Mee current Unit Activity 

This model considers four unit activities: ENGAGED, 

DISENGAGED, chemical decontamination (DECON) and 
resupply/refit (REFIT). ENGAGED describes a unit involved in 
a generic direct/indirect fire battle. DISENGAGED describes 
a unit that is in an assembly area or in transit. A unit in 
DECON is being cleaned of chemical contamination. A unit that 
is being resupplied or having maintenance performed on its 


equipment is in REFIT. 


ay 


2. ~~ 6Clurrent Unit Seeengen 

This model considers two factors that determine a 
unit's ability to fight and survive. The first factor is the 
amount of personnel, equipment and supplies the unit has on 
hand. This factor is a measure of the unit's MTOE (Modified 
Table of Organization and Equipment) strength, and is called 
a unit's MTOE factor. It takes on values of HIGH (minimal 
attrition losses to personnel and equipment), MEDIUM (moderate 
attrition losses to personnel and equipment) and LOW (severe 
attrition losses to personnel and equipment). The second 
factor further describes the condition of the unit's personnel 
and is called the PERSEFFECT (personnel effectiveness) factor. 
PERSEFFECT refers to the ability of the unit's personnel to 
perform their assigned tasks and is a function of the unit's 
exposure to chemical attacks. This factor has the values of 
HIGH (no degradation, personnel are totally effective) and LOW 
(personnel have been seriously degraded by exposure to 
chemical attacks). Determination of these categorical values 
is discussed in detail in section E. 

3. Current Unit Environment 

This refers specifically to the presence or absence 
of chemical weapons and protective gear. The three values for 
this indicator are CONVENTIONAL, CHEMICAL-MOPP (Mission 
Oriented Protective Posture) and CHEMICAL-NO MOPP. 


CONVENTIONAL refers to a state in which the unit is not 
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currently exposed to chemical hazards. A unit in CHEMICAL- 
MOPP is exposed to chemical hazards but has the required 
chemical protective gear to forestall loss of life strictly 
due to chemicals. CHEMICAL-NO MOPP refers to a unit exposed 
to chemical hazards that either does not have chemical 
protective gear or whose personnel have not yet donned their 


protective gear. 


D. MODEL STATE SPACE 
1. Number and Type of States 

Using the terms defined in section C. the model 
consists of 48 states numbered (0-47} with 46 transient states 
and two absorbing states. In actual combat units are rarely 
absorbed unless captured or destroyed by catastrophic losses. 
Replacements and reconstitution of weak units result in units 
surviving as entities even though all or most of the original 
personnel and equipment are gone. For the purposes of this 
model absorption means that the unit must go to the next 
higher level of restorative processes and is lost to the 
current level for a long enough time to be considered 
absorbed. The two absorbing states are COMBAT INEFFECTIVE- 
ATTRITION (unit has lost enough MTOE equipment/personnel to 
be totally combat ineffective) and COMBAT INEFFECTIVE- 
DEGRADATION (personnel are so degraded as to be totally 
ineffective). Appendix A contains the complete listing of the 


model state space. 
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2. Basic Model Assumptions 

Appendix B contains the model connectivity matrix that 
shows all the allowable one-step transitions. The basic 
assumptions used to determine what are allowable transitions 
are listed below with explanations. 

a. No Catastrophic Losses 

Units do not go from HIGH levels for either MTOE 

or PERSEFFECT directly to absorbing states. Units will 
transition through a series of states where these strength 
indicators are lowered one level at a time. Transitions 
directly to absorbing states are allowed only from LOW levels 
of either strength indicator. A model improvement would be 
to use actual random sojourn times and allow for a positive 
probability of very short or instantaneous transition times. 
Such a model would allow for units to transition quickly 
(possibly instantaneously) from HIGH levels to INEFFECTIVE 
levels (absorption) thus capturing the effects of catastrophic 
losses. 

b. Decontamination Before Refit 

If a unit has contacted chemical hazards it must 

first be decontaminated before it is allowed to go to refit 
and get supplies and replacements for attrition losses. This 
is in keeping with current doctrine so that the supply and 
maintenance facilities will remain uncontaminated and 


efftetence 


Zao 


c. MTOE Losses 
This model does not consider the personnel and 
equipment losses that occur through disease, accident and 
wearout while a unit is not actively engaged with the enemy. 
However, the MTOE factor of a unit can be lowered by any of 
the CHEMICAL-NO MOPP states because of the possibility of 
chemical casualties to personnel not in chemical protective 
gear. 
ad. Degradation While Disengaged 
In contrast to MTOE losses, this model does 
consider the loss of personnel effectiveness while units are 
not in active combat. Personnel will continue to suffer 
degradation as long as they are exposed to chemical hazards. 
The rate at which personnel become degraded is faster for 
ENGAGED units because of the greater exertion required in 
battle. Also, once they have been decontaminated, personnel 
will continue to be degraded unless they are given time to 
rest and recover. This rest and recuperation is handled in 
the model by having longer DECON sojourn times for units that 
Gomirom a LOWePERSEBFECT to a HIGH PERSEFFECT factor than for 
units that are cleaned up in DECON but maintain a LOW 
PERSEFFECT indicator. 
e. MOPP is 100% Protective 
Personnel exposed to chemical hazards while 
wearing MOPP will not be assessed casualties but the 


PERSEFFECT indicator may be lowered on the next transition. 
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f. Disengage Before/After Restoration 
Units must be in a DISENGAGED state before they 
can transition into the restorative processes, DECON or REFIT. 
Also when units leave DECON or REFIT they will go to 
DISENGAGED states. 
3. Interpretation of Transition Probabilities and Times 
The individual transition probabilities, p,;'S, can be 
considered a mix of decisions and externally imposed events. 
For example, if given a choice, a unit would have a higher 
probability of going to a disengaged state from an engaged 
state with both LOW strength levels, than from an engaged 
state at both HIGH strength levels. Thus the transition 
probabilities in this case reflect a conscious choice. (If the 
tactical scenario will allow a choice). On the other hand, 
the probability that a unit will receive a chemical attack as 
opposed to a strictly conventional weapon attack is dictated 
by the scenario and can be considered an external event 
imposed on the unit. 


The individual transition times, t.,'s, also have dual 


iJ 
interpretations. The first is that the transition time is 
the time it takes for the event that causes the unit to 
transition from state 7 to state j. The transition times can 
also be varied by the experimenter to represent a capability 


or efficiency. For example, allowing a unit to remain engaged 


longer and still be at HIGH MTOE strength is an increased 
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capability while taking longer for a REFIT or DECON operation 
denotes a less efficient REFIT or DECON process. 

For the purposes of this thesis the transition 
probabilities and times were determined by the author based 
on military judgement and the basic model assumptions 


Guecussea in®section D. 2. 


E. CATEGORICAL ASSESSMENTS 

One of the key issues in this model is the assigning of 
HIGH, MEDIUM and LOW strengths for the MTOE and PERSEFFECT 
factors of a unit. What exactly are these levels and how do 
they affect the probability of a unit transitioning into other 
states? 

A thesis by Cpt. Paul Crawford, Dynamic Study of Factors 
Impacting on Combat Power, addresses this issue. In his 
thesis Crawford considers a unit whose state is described by 
the percentage of personnel(PER), ammunition(AMMO), combat 
vehicles(VEH) and POL( Petroleum, Oil and Lubricants) on hand 
as opposed to the authorized unit levels for these variables. 
His problem was to map this state vector into a measure called 
the "unit effectiveness". The unit effectiveness has five 
levels: Totally Effective, Effective, Marginally Effective, 
Ineffective and Totally Ineffective. Personnel and combat 
vehicles were evaluated at three levels (100%, 75%, 50%) and 


ammunition and POL were evaluated at four levels (100%, 75%, 
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50%, 25%). Throughout the experiment the unit mission and 
threat were held constant. [ Ref. 4:p. 19] 

The experiment consisted of a questionnaire given to 45 
Army officers who then made a categorical judgement as to the 
unit's effectiveness level for a given state vector. The 
categorical responses were then transformed to an interval 
scale and the different state vectors were mapped into the 
interval scale. Crawford determined that the final function 
that mapped a state vector into the interval scale and thus 
into a category of effectiveness was elliptical of the form 


given by (7): 


Y= a° + b°(PER - X)? + c? (AMMO - X)? + a? (VEH - X)* + ©? (POL - X)*. (7) 


The best fitting model as determined by Crawford was: 


Y = 88.978 - .0056xX] -— .0055xX2 - .0054xX3 - .0005xX4, (8) 


Y = The value of unit effectiveness [interval scale], 


¥1 ="(Per — aoa) -: 
X2 = (AMMO - 100)‘, 
X3 = (VEH - 100)° ana 


X4 = (POL - 100)°. [ Ref. 4:pp. 11-23] 
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In the semi-Markov model considered in this thesis the 
state variable becomes the effectiveness (categorical) for the 
levels of HIGH, MEDIUM, LOW and INEFFECTIVE. The MTOE 
effectiveness in this model is exactly analogous to the unit 
effectiveness considered by Crawford. The PERSEFFECT 
effectiveness could also be calculated with Crawford's method 
but the state vector used would be different. The DEGRADATION 
state vector would include, type of chemical agent used, 
temperature, exertion level of troops, protective equipment 
used and time spent in the protective gear. 

This section has presented a method and reference for the 
user to determine what constitutes a HIGH, MEDIUM, LOW and 
INEFFECTIVE level for a unit. Even though categorical 
assessments are subjective, Crawford showed in his thesis that 
a group of decision makers can agree on a range of state 
vectors that map into a specified level of unit effectiveness. 
Also, by being able to fit a model to the results of the 
decision makers, Crawford showed a common systematic approach 
that was consciously or unconsciously used by them in their 
assessment. Thus if the method used for the categorical 
assessments, HIGH, MEDIUM, LOW and INEFFECTIVE, are explicitly 
stated then different users can use this model and assign 


subjective transition probabilities and times. 
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F. DEMONSTRATION OF MODEL 

In this section the working of the model and the output 
available from the model is demonstrated. A concern 
throughout the model development is, whether the model behaves 
reasonably and gives results that are believable. Assume that 


a unit can start in one of the following initial states: 


1. State 0: ENGAGED-CONVENTIONAL~HIGH-HIGH 

2 State 9: DISENGAGED-CONVENTIONAL-MEDIUM-LOW 
Se State 12: ENGAGED-CHEMICAL-~MOPP-HIGH-HIGH 

4. State 19: ENGAGED-CHEMICAL-NO MOPP-HIGH-LOW 
5. State 35: DISENGAGED-CHEMICAL-NO MOPP-LOW-LOW 
6. State 44; REF IT=-LOW=HIGH 


Figure 2 shows the expected time until absorption for the 
unit given the various initial states and transition 
probabilities and times given in Appendices C and D. 

The unit lifetimes vary from a high of 38.8 Hrs to a low 
of 13.8 Hrs. The lifetimes also follow a relatively intuitive 
pattern with the longest expected lifetime for the unit in 
REFIT and the shortest expected lifetime for the unit that is 
in chemical combat without protective gear and already at a 
LOW PERSEFFECT level. 

Reporting the expected lifetime is only part of the story. 
Also important is what was the unit capable of doing while 
"alive"? To gauge how a unit spent its useful life, the 


following classes of states were considered. 
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For the purposes of this thesis these are the classes of 
states considered. Now that the expected unit lifetimes have 
been determined the question is, what did the units actually 
do before they became combat ineffective? Figure 3 shows the 
percentage of time the unit spent in each of the above classes 


of states for each of the six initial states in this example. 


Unit fife Distributions for Six Initial 
states 
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Figure 3. Expected Distribution of Unit Lifetime for 
Various Initial States. 


Note that starting in state 44 gives the longest expected unit 
lifetime but that 40% of its life is spent in REFIT. This 
highlights the necessity to report expected unit lifetimes in 
conjunction with a breakout showing were the unit spent its 


useful life. Actually the times spent in COMCLEAN for the 
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initial states of 44, 0 and 6 are approximately equal, being 
between 5.7 to 5.2 hours for all three cases. 

The last area to address is the cause of the unit becoming 
combat ineffective. Did the unit become COMBAT INEFFECTIVE- 
ATTRITION (ATTRIT) or COMBAT INEFFECTIVE-DEGRADATION (DEGRAD) ? 
Using this model one cannot say what caused a specific unit 
to become combat ineffective but we can give a probability 
that a unit was absorbed in one or the other of the absorbing 
states. Figure 4 shows the probabilities of absorption for 


the six initial states in this example. 


Combat Ineffective due to Attrition vs 
Degradation for Six Initial States 
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As one would expect, the initial state of ENGAGED- 
CHEMICAL-NO MOPP-HIGH-LOW (State 19) has a very high 
probability of absorption in COMBAT INEFFECTIVE-DEGRADATION 
as opposed to COMBAT INEFFECTIVE-ATTRITION. 

An example of the output used to develop these results is 
contained in Appendix E. The summary report gives the 
information found in Figure 2 and Figure 3. A more detailed 
report is also provided by the computer program. The detailed 
report shown in Appendix E. gives the number of visits, 
expected sojourn time per visit and the total time spent in 
each state for which the unit spends at least one percent of 
its total life. The solution time required to run one case 
(a fixed P and T matrix) 1S approximately 3-5 minutes uSing an 


80286 based personal computer. 


G. MODEL USE AND CAVEATS 

Because this 1S an analytical model, the user needs to 
understand the abstractions and limits of the model. First, 
unlike a time step simulation, this model does not track an 
individual unit throughout its lifetime. For example a unit 
that "lives" nine hours could not possibly visit 23 states 
when travel times alone between activities are greater than 
nine hours. What the model is providing iS a nine hour 
expected lifetime for all units that start in the given 
initial state. Also, if a large number of units start in the 


Same initial state, we expect they will spend their lifetimes 
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distributed among the 23 states in the given proportions. 
Second, the expected time until absorption must always be 
considered in conjunction with the distribution of states in 
which the unit spent its useful life. Also, because the user 
makes the categorical judgments as to what constitutes a HIGH, 
MEDIUM and LOW strength unit any use of the model must 
consider how the judgments were made. Finally, because the 
user sets the P and T matrices the model is very scenario 


specific as are most Simulations with user input parameters. 


Sy IL 


V. ANALYSIS USING MODEL 


With the model in place, sensitivity analysis and testing 
of various assumptions can now be done. This analysis will 
concentrate on the three questions presented in the 
introduction : estimated threat, availability and priority of 
decontamination, and assumed degradation rates for personnel 
in a chemical environment. 

The following five step structure will be used to examine 
each of the questions. First, pose the question as a tactical 
Gecision maker would. Second, identify the MOES (Measures of 
Effectiveness) that answer the question. Third, describe what 
parameters in the model will be varied for the test. Fourth, 
present and analyze the raw model output. Fifth, frame the 
answer to the question in terms that are useful to a tactical 


decision maker. 


AS EST EMATE DO ine 
1. Question to be Answered 
The first question considered is, how does the 
frequency of chemical agent employment by the enemy affect the 
survivability of our force? The frequency of enemy chemical 
usage could be varied independently for each of two types of 
attacks: chemical weapons only and combined 


chemical/conventional weapons attacks. This analysis will 
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consider the relative frequency of all attacks that involve 
chemicals whether or not conventional weapons are involved. 
2. Supporting MOEs 

Several Measures of Effectiveness (MOE) address the 
somewhat vague term of survivability. The first MOE will be 
the time until a unit becomes combat ineffective given an 
initial starting state. The second MOE is the percentage of 
its useful life the unit spends in each of the six classes of 
states outlined in Chapter IV. The third MOE is the 
probability of absorption in each of the two absorbing states 
(i.e. the probability a unit becomes combat ineffective due 
to ATTRITION as opposed to DEGRADATION). 

For purposes of this analysis all units will start in 
state 6 DISENGAGED-CONVENTIONAL-HIGH-HIGH (owes, mot 
fighting, in an environment free from chemical hazards, at 
full MTOE strength and with individual personnel at full 
effectiveness (not degraded)). Using the first MOE, time to 
absorption (unit lifetime), by itself can give misleading 
results. What is important is not only how long a unit 
Survives but also what a unit was capable of doing while it 
was alive. A unit that survives 6 days but spends all its 
time in refit and decontamination is much less useful than a 
unit that survives 3 days but spends all its time in combat 
inflicting losses on the enemy. A unit's time to absorption 


must be reported in conjunction with the time and percentage 


oe 


of time a unit spent in the six classes of states from Chapter 
LWe 

Once a unit has lived out its useful life, what 
finally causes the unit to become combat ineffective; 
degradation or attrition losses? To the next higher level 
unit that needs to restore the combat ineffective unit this 
information is vital. Knowing what proportion of times units 
will become combat ineffective for each cause, the supporting 
unit(s) can determine the proper mix of restorative processes 
needed to support subordinate units. 

These three MOES will be used to evaluate the 
survivability and combat effectiveness of the unit throughout 
the analysis presented here. The primary MOEs are time to 
absorption and the proportion of time the unit spends in each 
of the classes outlined above. The secondary MOE is the 
probability of absorption in either of the two absorbing 
states in this process. 

3. Experimental Design 

To simulate various levels of chemical usage by the 
enemy the following procedure is used. All transition 
probabilities for transitions that lead from CONVENTIONAL 
states to CHEMICAL states are varied by a common factor from 
the base case. Using a common factor for all the transition 
probabilities keeps the relative frequency of each type of 


chemical attack constant while increasing or decreasing the 
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overall enemy usage of chemicals. The five variations tested 
and their relation to the base case are listed in Table II. 


For the present base case chemical weapons are used 


oo 


20% of the time if the unit is at MEDIUM-HIGH or better and 


oe 


10% of the time if the unit is at MEDIUM-LOW or worse. The 


worst case tested assumes a chemical usage rate of 40% of the 
time if the attacked unit is at MEDIUM-HIGH or better and 20% 


of the time if the attacked unit is at MEDIUM-LOW or worse. 


Table II TEST CASES FOR THE CHEMICAL THREAT ASSESSMENT 


Case Name Variation from Base Case 


Base None (see base case 
in Appendix C.) 


Nochem All transition probabilities 
from clean to dirty states 
are set to zero. 


Chemonly All transition probabilities 
from disengaged to engaged 
states are set to zero. 


Halfchem All transition probabilities 
from clean to dirty states 
are set to one half the 
base case value 


Twicchem All transition probabilities 
from clean to dirty states 
are set to twice the 
base case value. 


4. Analysis of Model Output 


The effect of various assumed chemical threats on the 


time until a unit becomes combat ineffective is shown in 
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Proure or Under a conventional weapons-only scenario the 
expected life of a unit is approximately 69 hours. In a 
chemical weapons-only scenario the expected life of a unit is 
reduced to approximately 43 hours. The minimum life expectancy 
of a unit, 26 hours, occurs at a chemical weapons usage which 
is twice the current base case. The worst case expected 
lifetimes are 38% of the conventional weapons-only and 59% of 


the chemical weapons-only lifetimes. 
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Figure 5 Expected lifetime of a unit under assumed chemical 
usage rates. 





The worst case scenario assumes a chemical usage rate 
of 40% if the attacked unit is at MEDIUM-HIGH or better and 
20% if the attacked unit is at MEDIUM-LOW or worse. The major 


consideration seems to be the introduction of chemical weapons 
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into the battle. Once they are being used, a variation of 
400% in their frequency of use results in only a 38% variation 
in the expected lifetime of a unit. Figure 6 shows the 
distribution of its lifetime a unit spends in the classes of 


states mentioned in Chapter IV. 


Distribution of unit life for the 


different chemical threats. 
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Figure 6 Distribution of unit life for chemical threat 
cases. 





Finally, what causes the unit to become combat 
ineffective, attrition losses or degradation due to the use 
of chemicals and chemical protective gear? Figure 7 is a 
comparison of the probabilities of absorption due to attrition 


and degradation. Note there is a positive probability that 


3/ 


units will become combat ineffective due to attrition losses 
in the chemical weapons-only case. This is the result of 
chemical attacks on unwarned, unprotected units that result 


in both attrition loses and degradation. 


Probability of combat ineffective due to 


attrition versus degradation. 
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Figure 7 Absorption probabilities for the various chemical 
threat cases. 





5. Discussion 
To achieve the minimum expected lifetimes for friendly 
units, an enemy's best tactic is to use a combination of 
chemical and conventional weapons. Using conventional or 


chemical weapons exclusively does not result in the shortest 
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expected lives for friendly units. This result is not 
unexpected because in U.S. and Warsaw Pact doctrine, chemical 
weapons are seen as combat multipliers, weapons that may have 
a synergistic effect in conjunction with other weapons. This 
is accounted for in the model by having the fastest 
transitions to worse states for units involved in combined 
chemical/conventional battles. What is interesting is the 
appearance of an apparent optimal mix of chemical and 
conventional weapons. To this author's knowledge, no studies 
have been done using large simulation models to determine the 
existence or actual value of such a chemical to conventional 


weapons mix. 


Be PRIORITY AND AVAILABILITY OF DECONTAMINATION 
1. Question to be Answered 
This is a two part guestion with both parts affecting 
the survivability of the force in this model. First, how does 
the decision to let a unit fight dirty or go immediately to 
decontamination affect the life time of the unit and where 
the unit spends its useful life? Second, given a base case 
of decontamination decisions, how does the lack of 
decontamination assets affect the life of a unit? 
2. Supporting MOEs 
The MOEs are the same ones used to answer the previous 
question about the frequency of enemy chemical agent use: time 


to absorption, portion of time the unit spends in each class 


iN, 


of states and probability of becoming combat ineffective due 
to attrition losses as opposed to chemical degradation. 
3. Experimental Design 

The second part of this question is the more 
straightforward one to answer. Starting with the base case 
probabilities of entering decon, decrease all these decon 
probabilities by a common factor until no decon assets are 
available. The various scenarios for this analysis are listed 
in Table III. 


Table III SCENARIOS FOR THE AVAILABILITY OF DECON ASSETS 


Case Variation from Base Case 
Base None. (see Appendix  C.) 
Bais All probabilities 


decreased to 75% 
of base case. 


D50% All probabilities 
decreased to 50% 
of base case. 

D25% All probabilities 
decreased to 25% 
of base case. 

Dnone All probabilities 


decreased to 0O. 
No decon available. 


The first portion of this question concerns when the 
decision to decon a unit is made based on the unit's strength 


and degradation level. A unit must be disengaged in order to 
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proceed to decon and in actual combat may be restrained from 
going to decon due to the tactical situation. For the 
purposes of this analysis the probability of going to decon 
from various disengaged states is considered as a decision 
variable. The various decon decision scenarios are given in 
Table IV and Table V. 


Table IV DECON DECISION SCENARIOS 


Case Name and States Transition 
Description to,from Probability 
old new 


BASE (see Appendix C.) 


DLOW 24,36 4 0 
(only decon 26,38 4 0 
units that 28,40 nS 0 
are LOW 30,36 23 0 
degradation) S250 see 0 

34,40 AS: 0 


(all other probabilities remain the same 
as BASE case) 


DNMP 24,36 ~4 oe, 
(shift prlorieyeeco, 5 / 4 ne 
of decon to 26,38 4 a2 
units w/o ay, 39 6 73 
chem protective 28,40 ale - OS 
gear). 29,41 any MoS 
30,36 .3 6 
37 is 6 
32,38 2 4 
See Soe, ao 6 
34,40 e 6 
35,41 58 6 


ee aie 
4. Analysis of Model Output 


The supporting MOEsS are the same as for the threat 
analysis and the model results for the various decontamination 


availability cases are given in Figures 8, 9 and 10. As decon 
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Table V FURTHER DECISION SCENARIOS 


Case Name States Transition 
to, from Probabilities 
old new 
DHIGH 24,36 4 oi 
(conserve Co, 4 6 
fighting abidlmerye2i6728 ~4 4 
by shifting 27,39 =c a6 
decon priority 28,40 a5 eS 
to higher 29,41 nee or 
strength units). 30,36 ag 24S 
Spleen: 3 i, <3 245 
32°38 Ae 52 
cic Ect:) 3 3 
34,40 3 215 
35,41 3 215 
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becomes less available the units experience shorter expected 


lives. Also, the decrease in Unit lifetime until comoat ineffective 
for various levels of decon availability 


unit lifetimes is’ nearly pe 
linear with a loss”7 of | 


approximately five hours for i" 


every 25% decrease in the 
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unavailable at this level the probability that a unit goes 


ineffective due to degradation approaches 0.62 (Figure 10). 
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Figure 10 Unit combat ineffective due to attrition versus 
degradation for various decon decision scenarios. 





For the decision cases the model output is shown in 
Figures 11, 12 and 13. The model shows the unit lifetimes and 
distributions to be very robust against changes in the decon 
priority, with approximately a 27% difference between the high 
and low expected unit lifetimes. The optimum scenario 
(longest expected unit lifetime) is the BASE (original 


decontamination decision) case. 
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1SIiOn Scenarios. 


Figure 12 Unit life distribution for various decontamination 
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In the current model the units are getting a double penalty 
because of the lack of decontamination assets. An underlying 
assumption is that units must be clean to go to refit and get 
supplies and replacements for attrition losses. 
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Figure 13 Absorption probabilities leone various 
decontamination decision scenarios. 


This is in keeping with current doctrine so that the supply 
and maintenance facilities will remain clean and efficient. 
A commander faced with no refit at all or doing it "dirty" 
will undoubtedly designate a lucky unit to perform the dirty 
refits and accept the Pon amneeon dimg=y BUMEcS iing in. 
This can be reflected in the model by the creation of new 


dirty refit states that have a longer turnaround time. Also 
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units undergoing dirty refit would have a small probability 
of becoming ineffective due to degradation. The restoring 
unit would also have a positive probability of going 


ineffective due to degradation. 


C. MOPP GEAR DEGRADATION EFFECTS 
1. Question to be Answered 
How does the assumed degradation that MOPP gear 
imposes on the personnel in a unit affect the unit's 
performance? 
2. Supporting MOEs 
The time until a unit becomes combat ineffective and 
the distribution of states in which the unit spends its useful 
life are the MOES for this analysis. Probability of 
absorption will not be considered for this question. In the 
experimental design the transition probability matrix is not 
changed, therefore the probabilities of absorption will remain 
constant for all cases considered. 
3. Experimental Design 
The degradation effects of a chemical environment are 
assumed to cause a loss in the unit's ability to sustain 
attrition and degradation strength levels as compared to a 


non-chemical environment. Both the degradation caused by the 
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chemicals themselves and the degradation due to the wearing 
of MOPP gear are considered.’ 

In this model these decreased capabilities are 
represented by decreased transition times to the next lower 
strength level. For example a unit that is HIGH-HIGH and in 
conventional combat can sustain HIGH-HIGH for one hour before 
going to MEDIUM-HIGH. Also this unit can disengage at 75% of 
it's HIGH to MEDIUM transition time and maintain it's HIGH 
MTOE strength rating. In contrast, a unit that is engaged and 
in MOPP will transition to MEDIUM in one hour but must 
disengage at 60% of this transition time to maintain it's HIGH 
MTOE rating. Also, if the PERSEFFECT indicator is already at 
LOW the unit's transition time to MEDIUM MTOE is .6 hours. 
For units that are engaged in chemicals and are not wearing 
MOPP the transition times are one half the MOPP transition 
times. 

This test used the BASE case transition time matrix 
(the TBASE case) and two other boundary cases, TBAD and TGOOD. 
In the base case units in CONVENTIONAL states can sustain 
strengths longer than units in CHEMICAL-MOPP, which in turn 
can sustain strength levels longer than units in CHEMICAL-NO 
MOPP. 


Degradation effects of chemical are those effects that 


hamper personnel but do not require medical attention to 
preserve the life of the soldier. Such effects include 
pinpointing of the pupils (miosis, nerve agents) and minor 
blistering from mustard agents. MOPP gear degrades troops by 
sensory deprivation and thermal stress (these suits get HOT). 
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In the TBAD case the transition times for CHEMICAL- 
MOPP units are 1.1 times the transition times for CHEMICAL-NO 
MOPP units. This represents the case were MOPP gear either 
does not protect adequately or is so stressful for troops to 
wear that it is only slightly better than no _ chemical 
protective gear. 

In the TGOOD case the transition times for CHEMICAL- 
MOPP units are 90% of the transition times for CONVENTIONAL 
Wnts. This represents the case were MOPP gear protects 
adequately and only slightly restricts troops in their mission 
performance. Also for TBAD and TGOOD, the transitions while 
DISENGAGED that have no counterpart in the CONVENTIONAL states 
are decreased and increased appropriately. 

4. Analysis of Model Output 

The assumed degradation rates were tested against the 
aifferent chemical threat levels described in Table II and the 
results are shown in Figures 14, 15 and 16. Figure 14 shows 
that as the assumed degradation rate increases (transition 
times decrease) the expected unit lifetime decreases. Though 
the lifetimes are affected more in the CHEMONLY case, the 
model does not exhibit any clear breakpoints, a case where the 
change in assumed degradation rates makes a dramatic or 
unexpected change in the expected unit lifetime. The effects 
of the assumed degradation rates on the unit lifetime 


distribution are shown in Figures 15 and 16. 
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Figure 14 Unit expected lifetimes for various assumed MOPP 
degradation rates. 





The degradation rates affect the CHEMONLY case much 
more than they affect the BASE threat case. Note that under 
TBAD, highest degradation effects, a unit spends 53% of its 
time in decontamination (see Figure 16). For the CHEMONLY 
case and TGOOD time matrix the unit spends only 41% of its 


time in DECON and the time spent in DISDIRTY states jumps from 


at 


28% to 45%. Note that for the BASE case the distribution of 
time spent in the different classes of states is almost 
constant with the only change in ranking occurring between the 


DISDIRTY and DECON states (see Figure 15). 
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Figure 15 Comparison of unit life distribution for the BASE 
threat case and three different degradation rates. 


5. Discussion 
In this case the model provides consistent results. 
If personnel are degraded worse the unit can expect to become 


combat ineffective quicker and this effect is more pronounced 
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for higher rates of enemy chemical weapons usage. Also, if 
MOPP gear is more efficient troops can fight in it longer and 
spend less percentage of their time in DECON. The model also 
is very robust to changes in the transition times. Some of 
the transition times are varied by a factor of five yet the 
expected unit lifetime and life distributions do not exhibit 


jumps of that magnitude (see Figures 14, 15 and 16). 


Unit. Hite distribution for @6HEMONLY 
threat and 3 assumed degradation rates. 


MB cOomMcLEAN 
KY compI ARTY 
([M} pi scLean 
B93 DISDIATY 


lite 
0:6: oe. “4 
stetstetatens 


OL 
( 


+ 6.0, 0, 640,¢ 
RRR ooo 
6.600 6:0 
£25525 50505 


* 

6 
Ox 

D505 


6,%,.¢ 
BESSA 
wart atatet 
eens 
Motete% 
€:0:0:6:0 
eeecee 


© 
Oe 


eee 
(yy 
6 
®& 

-0-0,0,0° 

2250505 


© 
9%, 
O 
, 
O 


‘Cs 

(\ 

e 
6:0:6°6 


+ 
O 
O 


é 
4 
6 
o.! 
550525 


. 

POO 

OOO) 
? 


~ 
Cc 
> 
— 
©) 
yw, 
e 
@ 
O 
§ 


2 
\/ 

POARKR 
405 


e@ 
6 
524 


O 
»~) 
eo 
ca 29% 2 
nheteceansn a ascaeet 
6 
005% 2, 


is 
a 


© 


4 
e “6 
Ox) 


aa 
ore 
© 


stetete' ROP 


® 
©) 
one: 
Sy, 
m6 
& 
at yte' 


oO 

e& 

rim 
© 


ee 

O,8 

©) 
aM a%a®, 


Ss 


O 
@e@ 


= eS, 
OM) 
rere” 
slettate 

er etek 


oe 


O 
2 


O 
CRM R RR R 


O 


0;0.9, 
Q 
OX) 
We 


O 


0,98 @. 
QS 
\/ 
& 
oor 
os 
*e; 
2. 


CY 
OO 


s me ) 
$0 
ete! 


o 
Para’ 
| VaVat, OC, CpG aH, 
e058, 
©) 
ras 


e 
e 


“oe 
= 
seco 


® 
e@ 
we 

2%, 
‘7 





Figure 16. Comparison of unit life distribution for the 
CHEMONLY threat case and three different degradation rates. 


D. SENSITIVITY ANALYSIS 
Sensitivity analysis is concerned with how general are the 


results from an experiment? Are the results robust, 
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applicable over a range of input parameters, or are the 
results so sensitive that one minor change in the inputs will 
cause entirely different conclusions? For this particular 
model the inputs are the P and T matrices. Because the 
experiment with the assumed degradation rates showed the 
expected unit lifetime to be less sensitive to changes in the 
T matrix than in the P matrix, a sensitivity analysis was 
performed on the P matrix (transition probabilities). 

The original threat-cases experiment was performed by 
systematically varying part of the P matrix for the model and 
holding the rest of the transition probabilities fixed 
relative to one another (see Appendix F). How sensitive are 
results of the previous analysis to changes in the portion of 
the model that are not being varied in a systematic manner? 
To test this sensitivity the following experiment was 
conducted. The five chemical threat levels were tested 
against the base P matrix and against five matrices that had 
been perturbed using a Normal(0,.01) random variable. The 
N(0O,.01) random variable was added to a matrix of all ones the 
same size as the model P matrix. The P matrix, with the test 
entries removed, was the multiplied by this randomized matrix. 
The resulting probabilities were then normalized so that the 
rows still summed to one when the removed probabilities were 
added back in. The results for the total unit lifetimes are 


shown in Figure 17. The randomized matrices caused no change 
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in the relative ranking of the threats but did affect the 
absolute lifetime of the unit particularly at the extreme 
cases of no chemicals and all chemicals. P10T1 is the 


Original transition probability matrix given in Appendix Cc. 
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Figure 17 Total unit lifetime for the base case and five 
randomized P matrices. 





Figure 18 and Figure 19 show the unit lifetime 
distribution for the NOCHEM and CHEMONLY threat cases. These 
were the cases whose expected unit lifetimes showed the most 
variance. The distribution of unit lifetime appears to be 


more robust than the expected unit life. 
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Figure 19 Distribution of unit lifetime for the chemical 
only case. 


57 


VI. CONCLUSION AND RECOMMENDATIONS 


The semi-~Markov model presented in this thesis has shown 
itself to be a useful tool for sensitivity analysis and the 
ranking of alternative threats and courses of action. The 
model provides the following insights into the three questions 
posed in the Introduction to the thesis. 

The major factor that determines the expected lifetime and 
capabilities of a unit is whether or not chemical weapons have 
been introduced into the battle. Units can expect to live 
about 70 hours ina strictly conventional battle whereas their 
expected lifetime drops to 40 hours with the introduction of 
chemical weapons at the lowest frequency of usage tested in 
this model. Once chemical weapons have been introduced, 
however, increasing their use does not continue to lower unit 
lifetimes as dramatically. A 400% increase in the use of 
chemical weapons results in approximately a 38% decrease in 
the expected life of a unit. As chemical usage increases, 
however, more of a unit's lifetime is spent in restorative 
processes and not in productive combat. Also, the model 
points to the existence of an optimal conventional/chemical 
weapons mix. 

As regards the availability of decontamination assets and 


their use, the model shows that the availability of assets 
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affects the survivability of units more than the decision of 
when to decontaminate a unit. The decontamination priority 
scheme chosen for the base case of the model gave the best 
results for expected unit lifetime by a small margin over the 
decision to decontaminate units at HIGH MTOE strengths. As 
might be expected, the more that chemical gear degrades 
individual soldiers, the shorter the expected unit life will 
be. The results of the degradation experiment show how robust 
the model is to changes in an individual transition time or 
even a group of transition times. Varying individual 
transition times by a factor of five did not vary unit 
lifetimes by even a factor of one-half. This robustness also 
points out that varying the transition times alone may not be 
an adequate experiment. Further experiments should consider 
varying both the transition times and transition probabilities 
to test the effects of assumed degradation rates. 

The model is not intended as a replacement for time step, 
high-resolution simulations. However, using such simulations 
as input for the P and T matrices the model is a quick way to 
answer " What if...?" questions concerning changes in the 
threat and changes in tactics. 

This class of analytical models is limited only by the 


user's patience in designing the state space and experiments. 
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APPENDIX A. STATE SPACE LIST 


The state space used for this model is given below. The 


terms are defined in Chapter IV Model Formulation. The bold 


numbers are the state number and the written portion is the 


state description. 


WOON AU PWN Ee O 


ENGAGED-CONVENTIONAL-HIGH-HIGH 
ENGAGED-CONVENTIONAL-HIGH-LOW 
ENGAGED-CONVENTIONAL-MEDIUM=-HIGH 
ENGAGED-CONVENTIONAL-~MEDIUM-LOW 
ENGAGED-CONVENTIONAL-LOW-HIGH 
ENGAGED-CONVENT IONAL=- LOW- LOW 
DISENGAGED-CONVENTIONAL~HIGH-HIGH 
DISENGAGED-CONVENTIONAL-HIGH-LOW 
DISENGAGED-CONVENTIONAL-MEDIUM-HIGH 
DISENGAGED-CONVENTIONAL-MEDIUM-LOW 
DISENGAGED-CONVENTIONAL-LOW-HIGH 
DISENGAGED-CONVENTIONAL=-LOW- LOW 
ENGAGED-CHEMICAL-MOPP-HIGH-HIGH 
ENGAGED-CHEMICAL-MOPP-HIGH=-LOW 
ENGAGED~CHEMICAL“~MOPP-MEDIUM-HIGH 
ENGAGED-CHEMICAL-MOPP-MEDIUM-LOW 
ENGAGED-CHEMICAL-MOPP-LOW-HIGH 
ENGAGED-CHEMICAL-MOPP-LOW-LOW 
ENGAGED-CHEMICAL=-NO MOPP-HIGH-HIGH 
ENGAGED-CHEMICAL-NO MOPP-HIGH-LOW 
ENGAGED-CHEMICAL-NO MOPP-MEDIUM-HIGH 
ENGAGED-CHEMICAL-NO MOPP-MEDIUM=-LOW 
ENGAGED-CHEMICAL-NO MOPP-LOW-HIGH 
ENGAGED-CHEMICAL-NO MOPP-LOW-LOW 
DISENGAGED-CHEMICAL-MOPP-HIGH-HIGH 
DISENGAGED-CHEMICAL-MOPP-HIGH-LOW 
DISENGAGED-CHEMICAL-MOPP-MEDIUM-HIGH 
DISENGAGED-CHEMICAL-MOPP-MEDIUM=-LOW 
DISENGAGED-CHEMICAL-MOPP-LOW-HIGH 
DISENGAGED-CHEMICAL-MOPP-LOW-LOW 
DISENGAGED-CHEMICAL-NO MOPP-HIGH-HIGH 
DISENGAGED-CHEMICAL-NO MOPP-HIGH-LOW 
DISENGAGED-CHEMICAL-NO MOPP-MEDIUM-HIGH 
DISENGAGED-CHEMICAL-NO MOPP-MEDIUM-LOW 
DISENGAGED-CHEMICAL-NO MOPP-LOW-HIGH 
DISENGAGED-CHEMICAL-NO MOPP-LOW-LOW 
DECONTAMINATION-HIGH-HIGH 
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DECONTAMINATION-HIGH-LOW 
DECONTAMINATION-MEDIUM-HIGH 
DECONTAMINATION-MEDIUM-LOW 
DECONTAMINATION-LOW-HIGH 
DECONTAMINATION-LOW-LOW 
REFIT-MEDIUM-HIGH 
REFIT-MEDIUM-LOW 
REFIT-LOW-HIGH 
REFIT-LOW-LOW 

COMBAT INEFFECTIVE-ATTRITION 
COMBAT INEFFECTIVE-DEGRADATION 


of 


APPENDIX B. CONNECTIVITY MATRIX 


The connectivity matrix for the model is given on the next 
page. The bold numbers are the state numbers and the entries 
in the cells denote whether or not a transition from state 7 
to state J is allowed in one transition. An "X" denotes a 
positive probability of transitioning from state 7, the row 
number, to state j, the column number. A blank denotes that 


the transition is not allowed or at least not considered in 


the present model. 
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APPENDIX C. ONE-STEP TRANSITION 
PROBABILITY MATRIX 
This appendix contains the the numerical values for 
transition probabilities that correspond to the "X"s on 
connectivity matrix in Appendix B. The probabilities 
broken down by rows and each section is labeled with 
current state, row, number. The notation p(N,M) is 


probability of going from state N to state M in 


the 


the 


are 


the 


the 


One 


transition. All probabilities not expressly assigned in this 


appendix are equal to zero for this model. 
State 0: p(0,2)= 0.2.6, p(0,6)= 0227 pie l2) — Ono 
p(0,18)= 0.05 


State 1: p(1;3)= 0.5) pil,7)=]= O23) pil ts) — soe 
jog UP PO), 2) 


State 2: p(2,4)= 0.3, pl2,8)= 0.5, pe, 4) =7 
p(2,20)= 0.05 


State 3: P(3,5)= 0.2, pl(3,9)= 0.77 pls) oor So 
bits, 21) = 20 70> 


State 4: p(4,10)= 0.7, p(4,16)= 07075, pia, 22 —) 6. 0257 
p(4,46)= 0.2 

State 5: p(5,11)= 0.6, p(5,217)=.0. 075 eG 2) eee 
p(5,46)= 0.3 

State 6: p(6,0)= 0.8, p(6,12)— 0.2075 pie io) ee On O72 oF 


p(6,24)= 0.075, p(6,30)= 0.025 


Statenv: P(7;1)="0.8, pl(7,13)— CeOy SC lc) pe ise 
D(7,25)= 0.075, p(7 oon > 


State 8: p(8,2)= 0.6, p(8,14 )="0.075, 3 b( a 20) =~ 30. e225. 
P(8,26)= 0.075," p(s; 32) —s0. ez. 
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State 9: Dice) — Ove roma — 0.0975, pro 2l)= 0.0125, 
p(9,27)= 0.0375), p(9,33)= 0.0125, p(9,43)= 0.3 


State 10: p(10,4)= 0.3, p(10,16)= 0.0375, p(10,22)= 0.0125, 
p(10,28)= 0.0375, p(10,34)= 0.0125, p(10,44)= 0.6 


State 11: p(11,5)= 0.3, p(11,17)= 0.0375, p(11,23)= 0.0125, 
p(11,29)= 0.0375, p(11,35)= 0.0125, p(11,45)= 0.6 


State 12: Dero) Oe ep i2y 14 0.2, pil2,.15)= 0.1, 
0 


p(12,24)= 0.6 
Siate. a. Digkoeel > oes Gis 25)— Osea, pri3,47)— 0.1 
State 14: p47 15)— 10 .2,) pcl4, 16)= 0-2, p(i4,17)= 0.2 
Dildo) — 10es 
State 15: Dio le) — Omeemo(ts,2/)— 0.6,6p15,47)= 0.1 


State 16: p(16,17)= 0.2, p(16,28)= 0.4, p(16,46)= 0.4 
Setee iy) p(i7,29)= 0.3, pll7,46)= 0.5, p(17,47)= 0.2 


State 18; p(18,12)= 0. y= 70, Pls ,20)= 0.2 
Dile..21) = ©. 1, .pi1s ,30)= 0.4 


Beare 19: p(l9i3j= O22, p(l9,21)= 0.4, p(l9,31)= 0.2 
p(19,47)= 0 


State 20: eZ Oe) — Be (20,219) — FOn28 pi2O,22)= 0.3 
DiGZOw2s )=— 02.2), p(20,32)= 0214 

State 21: preg lS = 6n2, Piz! 24) —30.4>, p(21,33)= 0.2 

p(21,47)= 0.2 


Stare 22: Pigevwloya 0.2), Pi22,,25)—70.2;, 9p (22,34)=° 0.2 
p(22,46)= 0.4 


State 23; D(23,.7)= 
p(23,47)= 


2) peZs es iO 2m p23, 46)— 0.4 


eee me On 242) — "ONS, pi24,25)= 0.3, p(24,36)= 0.4 
a eeme D259) — 064, pi 25,o7)—= 0.4, p(25,47)= 0.2 
State 26: p(26,14)= 0.3, p(26,27)= 0.3, p(26,38)= 0.4 
peony ee 27 lod) — 0.2) D(27,80)= 0.6, p(27,47)= 0.2 
State 28:  p(28,16)= 0.2, p(28,30)= 0.3, p(28,40)= 0.5 
State 29: 2 i) — Oe mez oa) — OnLy, p(29,47)= 0.2 
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State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


Stace 


State 


State 


State 


SOM 


oe 


32:3 


Son 


34: 


3S. 


Chere 


ay a 


Bion 


39° 


40: 


41: 


42: 


Ase 


44: 


45: 


46: 


47: 


p(30,18)= 0.2, 
D(30; 33 )=s022, 
p (31,19) =60 23 
p(31L,37)= 7@23s" 
Di 32420) =50 
p (32:;, 35) =90s25 
D(33721)— 0-32 
p(397 59) = ga" 
p(3427 22) — 0. 
p(34,40)= 0.3, 
p(35,23)= 028; 
Dita ,46)— Onez, 


p(36,6)= 1.0 

Dio7  6o)——Oe4e 
P(36, 5) — 0. 5, 
p(39,8)= 0.2, 
DAO a Ome 
D(41,,10)= "Ones 
Di426)— eo 

p(43,7)= 1.0 

p(44,6)= 0.6, 
OAS. = OT. 
p(46,46)= 1.0 


p(47,47)= 1.0 


p(30.24)= 
p(30,36)= 


p(31,25)= 
p(31,47)= 


p(32,26)= 
p(32,38)= 


p(33,27)= 
p(33,47)= 


p(34,28)= 
p(34,46)= 


p(35,29)= 
p(35,47)= 


D(37, 7 ).= 108 


p(38,42)= 0.5 


p(39 ,9)= 20. 
p(40,44)= 


p(41,11)= 


p(44,8)= 0O. 


p(45,9)= 0O. 
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6 


3, 


0; 


, p(30,31)= 


nee (315,33)= 


, P(32,33)= 


3, p(33,35)= 


6 


, p(34,35)= 


pubioo, 41)= 


p(39,43)= 0. 


-3;, PCr, 45) = 


APPENDIX D. TRANSITION TIME MATRIX 


This appendix contains the the numerical values for the 
transition times that correspond to the "X"s on the 
connectivity matrix in Appendix B. The times are broken down 
by rows and each section is labeled with the current state 
(row) number. The notation t(N,M) is the transition time for 
going from state N to state M in one transition. Note that 
times are only given for transitions with a positive 
transition probability from Appendix C. Times are expressed 


ious Ss. 


Seace 0: Oy 2) — se Opt 0, oO) — 9 O.7>, £0, 12)— 0.1, 
G(0718) = O27 

Stace 1: Cig sta Oe ot (de 7 a Ones 7 t( 1 3 j=*0.1, 
CCl 29)= 0.4: 

State 2: epee) — 2 Om ta(2 8) — doo, te (2,14)= 0.1, 
G(2,20)= 0f1 

Srace. 3: fee(e3 o>) eo.) = lage (ou, 15 )— Ol, 
to, 2h) = OC 

Stace 4; ako) — NO. 5s ete ae to)— Ole t(4,22)= 0.1, 
G(4,46)— 0.7 

State 5: (Sa) =O oe eto 7) On dete (5 ,23)— 0.1, 
€(5,46)= 0.5 

State 6: eco, Os 1. Or meete, 2 y—ae. OF, C(6,168)= 1.0, 


(Ep acy WW, 1 Cee et ate 


State 7: 7h) OCI o =O, (7,19 = 1.0, 
(7-2 oO C71 =P 1.0 


State 8: t(8,2)= 1.0, t(8,14)= 1.0, t(8,20)= 1.0, 
t(8,26)= 1.0, t(8,32)= 2.0 
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2 lis 
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Ls 
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25:5 


20. 


Pe or 
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29. 


t(9,3)= 1.0, 
t(9, 27). =m 


t(10,4)= 1.0, 
t(10,28)= 1.0, 


t(9,15)= 1.0, 
t(9,33)= 1.0, 


t(9,21)= 1.0, 
t(9,43)= 2.0 


ti(10 46) = ope 10, 22), 
t(10,34)= 1.0, 


t(10,44)= 2.0 


t(1lyg5)He020, tidy) =O eee) ao 
0, t(11,35)= 1.0, 


t(11,29)= 1. 


t(12,13)= 
t(12,24)= 


t(13,15)= 


t(14,15)= 
t(14,26)= 


t(15,17)= 
t(16,17)= 
t(17,29)= 


t(18,12)= 
aor: eee 


elon ies = 
t(19,47)= 


t(20,14)= 
t(20,23)= 


t(21,15)= 
t(21,47)= 


t(22,16)= 
t(22,46)= 


t(23,17)= 
t(23,47)= 


t(24,12)= 
t(25,13)= 
t(26,14)= 
t(27,15)= 
t(28,16)= 


t(29,17)= 


@ 1S) 


.75, 


Oo 


2) ©1527) = 102 


RA (Ae 
5 
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.6, t(13,25)= 0.36, 


.0, t(14,16)= 2.0, 
2 
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.O, t(16,28)= 1.8, 
3, t€(17,46)= Of537 


.2, t(18,19)= 1.5, 
t(18,30)= 0.25 


<2, 7t(19220) — Orer 
72, (20,200 =e 
t(20,32)= 0.5 


22, (2123) = Omer 


s27 (22529 —2 or 


35)= 0.1, 


.9, t(25,37)= 2.0, 
9, £(26,27)= 4. 0F 


9) t(27, 39) =oceor 


9, £(297 41 eon 


CC lr 4 y= 250 


0, €(12,14)= 1.0, Fee, 15)— eee 
6 


t(13,47)= 0.5 


+ ( Jide 


t(15,47)= 0.5 
t(16,46)= 0.7 
t(17,47)= 0.5 


t(18,20)= 0.5, 


t(19,31)= 0.15, 


(20), 22)= "1. OF 


t(21,33)= 0.3, 


t(22,34)= 0.75, 


£(23,46)= 0.2, 


9, t(24,25)= 4. Omer 6) = 2.0 


t(25,47)= 2.0 
t(26,38)= 2.0 


t(27,47)= 2.0 


.9, t(28,30)= 4.0, t(28,40)= 2.0 


£(29,47)= 2.0 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


State 


Cle 


33 


Bas 


33°3 


Sa: 


3.5.5 


36: 


3275 


33.5 


si o0 


40: 


Ani 


42: 


alee 


44: 


Ac 


46: 


47: 


t(30,18)= 0.5, 
EMeONss)= Oso. 
t(31,19)= 0.5, 
t(31,37)= 0.5, 
t(32,20)= 0.5, 
t(s2s5)= 1.0, 
t(33,21)= 0.5, 
13 3 On = Oy 


t(34,22)= 0.5 
t(34,40)= 2.0 


£(35,23)= 0.5, 
t(35,46)= 0.5, 


e136 ;,6)= 4.0 

(37 ,)6)= 6aor 
t(38,8)= 4.0, 
t(39,8)= 6.0, 
Econo — 440, 
ce eile (0))) Sl yee 
t(42,6)= 5.0 

t(43,7)= 5.0 

t(44,6)= 8.0, 
t(45,7)= 8.0, 
t(46,46)= 1.0 


t(47,47)= 1.0 


eS) 
t(30,36)= 


t(31,25)= 
t(31,47)= 


t(32,26)= 
t(32,38)= 


t(33,27)= 
t(33,47)= 


t(34,28)= 
t(34,46)= 


t(35,29)= 
t(35,47)= 


t(37,7)= 4. 


t(38,42)= 4.0 


t(39,9)= 4. 
t(40,44)= 


t(41,11)= 


t(44,8)= 8. 


t(45,9)= 8. 


69 


e) 


0, 


4. 


20; 


atte, 21)= 


, t(31,33)= 


, t(32,33)= 


mt (33 335) = 


mt (34,35)= 


, t(35,41)= 


OOS), ee - 


€(41,45)= 


2.0, 


APPENDIX E. SAMPLE MODEL OUTPUT 


This appendix contains the model outputs that were used 
for the sample runs in Chapter IV. The first series of 
reports are the Summary reports that are given for each 
initial state. The report shows the total time until 
absorption, the actual time (hours) spent in each class of 
states and the percent of the unit's total life that was spent 
in each class of states. Below the Summary report are three 
numbers; the initial state number, the probability of 
absorption in COMBAT INEFFECTIVE-ATTRITION and the probability 
of absorption in COMBAT INEFFECTIVE-DEGRADATION. The first 
set of reports are created by the function AUTO and are 
displayed by the function REPORT (see Appendix G. APL 


FUNCTICGNS) . 
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COFREAUw I) 4)10e? ) 
FOR THE BASE P10T1 T10T1 CASE 
AND INITIAL STATE OF ENGAGED-CONVENTIONAL-HIGH-HIGH 
THE TOTAL UNIT LIFE UNTIL COMBAT INEFFECTIVE WAS 29.840 HRS 


OF ITS USEFUL LIFE THE UNIT SPENT IT IN THE FOLLOWING AREAS 


HOURS PERCENT OF DUPFE 

COMBAT CLEAN 5.740 19.24 
COMBAT DIRTY a. 230 10. 85 
DISENGAGED CLEAN 5. 189 NS) 
DISENGAGED DIRTY 3.949 3223 
DECONTAMINATION 4.339 14.54 
REFIT 7.384 24.75 

0 0.5866583374 0.4133416626 


A 


SGOPREAD 2, 4909359 


FOR THE BASE PlOTa TIOTI CASE 
AND INITIAL STATE OF DISENGAGED-CONVENTIONAL-MEDIUM-LOW 


THE TOTAL UNIT LIFE UNTIL COMBAT INEFFECTIVE WAS 29.092 HRS 


OE ITS USEFUL, LIFE THE UNIT SPENT IT IN THE FOLLOWING AREAS 


HOURS PERCENT ©F Shire 

COMBAT CLEAN De Jo8 16 35 
COMBAT DIRTY URS 758! Dace 
DISENGAGED CLEAN volo 24.74 
DISENGAGED DIRTY 2002 Jaze 
DECONTAMINATION 3.647 Peon 
REerit 8.704 29.92 

9 0.4534389037 0.5465610963 


PORSTHEA BASE PLOTi- T1011 CASE 
AND INITIAL STATE OF ENGAGED-CHEMICAL-MOPP-HIGH-HIGH 


ie LOTAG UNET EIPeE UNTIL COMBAT INEFFECTIVE WAS 28.304 HRS 


OF ITS USEFUL LIFE THE UNIT SPENT IT IN THE FOLLOWING AREAS 


HOURS PERCENT OF EIFE 
COMBAT CLEAN a 4-29 2 = 2 
COMBAT DIRTY 4.109 14.52 
DISENGAGED CLEAN 4.084 14.43 
DISENGAGED DIRTY Sapo ne 19.47 
DECONTAMINATION Sl) S| 20.33 
REFIT 5.416 19.14 
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0.460368813 


0.539631187 


(OFREAD 12 4)(19;::;3 
FOR THE BASE P10T1 T10T1 CASE 


AND INITIAL STATE OF ENGAGED-CHEMICAL-NO MOPP-HIGH-LOW 


THE TOTAL UNIT LIFE UNTIL COMBAT INEFFECTIVE SHES 


13.768 HRS 


OF ITS USEFUL LIFE THE UNIT SPENT I7 IN THE FORTOWING AREAS 


HOURS PERCENT OF LIFE 
COMBAT CLEAN 17.416 26a 
COMBAT DIRTY 1.589 11.54 
DISENGAGED CLEAN 2.174 15279 
DISENGAGED DIRTY 2.200 15.98 
DECONTAMINATION 3 10010 Pls he 
REFIT 3.050 el Us 
192° 0.2929669279 0.70705 s0721 


COFPREAD de40 035 53) 
FOR THE BASE PrOrie 7 TOT VersSe 
AND INITIAL STATE OF DISENGAGED-CHEMICAL-NO MOPP-LOW-LOW 
THE TOTAL UNIT LIFE UNTIL COMBAT INEFFECTIVE WAS 22.178 HRS 


OF ITS USEFUL LIFE THE UNIT SPENT I) THSTHE FPOLLewiiGearens 


HOURS PERCENT OF LIFE 

COMBAT CLEAN 2-510 Lie s2 
COMBAT DIRTY ; oon 4.48 
DISENGAGED CLEAN 3.2094 Gee 
DISENGAGED DIRTY 3.022 MSO 2 
DECONTAMINATION 4.855 2.200 
REET 7.204 32.48 

pe 0.5294819827 0.4705180173 


(OFREAD 1 4)(44;33) 
FOR THE BASE -PlOT?] Tl0ti CASE 
AND INITIAL. STATE OF REFIT —-LOwW-pien 
THE TOTAL UNIT LIFE UNTIL COMBAT INEFFECTIVE WAS 22 oe 4enkS 


OF ITS USEFUL LIFE THE UNIT SPENT If IN THE FOrLroviINcee ras 


HOURS PERCENT OF LIFE 

COMBAT CLEAN sia) beh) 13.28 
COMBAT DIRTY S220 8229 
DISENGAGED CLEAN Go252 15.84 
DISENGAGED DIRTY A135 10.65 
DECONTAMINATION 4025 eS TGs 
REP 15.650 40.31 

44 0.5787679768 0.421 2320232 


oZ 


. 


This next report is a sample of the Detailed report also 
created by the function AUTO. This report shows the states 
visited, the number of visits, the sojourn time for each visit 
and the total time spent in each state for all states in which 
the unit spent at least one percent of it's total life. The 


heading information is the same as that for the Summary 


report. 
OFREAD 1 5 
FOR THE BASE P10T1 T10T1 CASE 
AND INITIAL STATE OF ENGAGED-CONVENTIONAL-HIGH-HIGH 

THE TOTAL UNIT LIFE UNTIL COMBAT INEFFECTIVE WAS 29.840 HRS 

OF ITS USEFUL LIFE THE UNIT SPENT IT IN THE FOLLOWING STATES 

STATES VISITS TIME EACH TOTAL TIME 

VISIT IN STATE 
ENGAGED- CONVENTIONAL-HIGH-HIGH 1.979 7 7G Tes26 
ENGAGED-CONVENTIONAL-HIGH- LOW .516 .674 .348 
ENGAGED-CONVENTIONAL-MEDIUM-HIGH O82 370 2.646 
ENGAGED-CONVENTIONAL-MEDIUM- LOW .569 1.350 .768 
ENGAGED-CONVENTIONAL- LOW-HIGH o7 oe e521 394 
DISENGAGED-CONVENTIONAL-HIGH-HIGH digees 1.000 Lace 3 
DISENGAGED-CONVENTIONAL-HIGH- LOW -645 1. 000 .645 
DISENGAGED-CONVENTIONAL-MEDIUM-HIGH ta 1.200 1.489 
DISENGAGED-CONVENTIONAL-MEDIUM- LOW OL Pes 0i0 - 674 
DISENGAGED-CONVENTIONAL- LOW-HIGH «590 £. 600 -944 
ENGAGED-CHEMICAL-MOPP-HIGH-HIGH .551 hs sate . 556 
ENGAGE D-CHEMICAL-MOPP-MEDIUM-H1IGH «O12 1.900 Lee 
ENGAGED- CHEMICAL-MOPP-MEDIUM- LOW .444 .842 Sele 
ENGAGED-CHEMICAL-MOPP- LOW-HIGH 7308 1.600 .493 
DI SENGAGED-CHEMICAL-MOPP-HIGH-HIGH .448 Daa 0 1007 
DISENGAGED-CHEMICAL-MOPP-HIGH- LOW .407 1.560 .635 
DISENGAGED- CHEMICAL-MOPP-MEDIUM-HIGH eo) 20 . 660 
DISENGAGED~CHEMICAL-MOPP-MEDIUM- LOW 994 7 8G . 696 
DISENGAGED~-CHEMICAL-MOPP- LOW-HIGH P52 2.380 .362 
DISENGAGED-CHEMICAL-MOPP-LOW- LOW .176 1.890 O32 
DECONTAMINATION-HIGH-HIGH 205 4.000 820 
DECONTAMINATION-HIGH- LOW er 4.800 .834 
DECONTAMINATION-MEDIUM-HIGH .126 4.000 .503 
DECONTAMINATION-MEDIUM- LOW 2a 4.400 gt SK) 
DECONTAMINATION-LOW-HIGH .083 4.000 e332 
DECONTAMINATION~LOW-LOW ey, 5.400 A740 
REFIT-MEDIUM-HIGH oi 5.000 1.555 
REFIT-MEDIUM- LOW m262 5.000 1.409 
REFIT- LOW-HIGH .404 8.000 320) 
REFIT~-LOW- LOW .149 8.000 1.189 
FF 
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APPENDIX F. EXPERIMENTAL IDENTIFIER MATRICES AND DATA 


This appendix contains connectivity matrices that identify 
which probabilities and transition times were changed for the 
experiments in Chapter V. The notation is the same as that 
used in Appendix B with the exception of special letters which 
are explained when the appear. 

The first matrix on the following page shows the entries 
which were set to zero for the CHEMONLY case. These entries 
correspond to transitions from DISENGAGED to ENGAGED states 
(conventional weapons attacks). The letter "E" denotes 
conventional weapons attacks only and "S" denotes simultaneous 
conventional and chemical weapons attacks. 

The second matrix shows the transitions between 
CONVENTIONAL to CHEMICAL states which were varied for the 
HALFCHEM, TWICCHEM, and NOCHEM cases. The letter "C'" denotes 
chemical weapons only attacks and "S" denotes simultaneous 
chemical and conventional weapons attacks. How they are 


varied is described in Table II., Chapter V. 
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Below is the conventional weapons attacks identifier matrix. 


1 2 3 4 
012345678901234567890123456789012345678901234567 
0 X X x x 
1 x X x X 
2 x x x x 
3 X x x x 
4 X X x x 
5 x X x x 
6E S S x x 
7 €E S S x x 
8 E S S x x X 
9 E S S x x x 
10 E S S x x x 
11 E S S X x x 
12 x X X X 
13 x x x 
14 x XX x 
15 x x x 
16 x x X 
17 x X X 
18 x Xx xX X x 
19 x X x X 
20 x x X X x 
21 x x x x 
22 xX X x x 
23 x x x xX 
24 E x X 
2 E x x 
26 E x x 
o7 E x Xx 
28 E x x 
29 e x X 
30 E x Xx xX x 
31 E X x x x 
32 E x x Xx x 
33 E x x x x 
34 E x x x X 
35 E x x x xX 
36 x 
37 x xX 
38 x x 
39 xX xX x 
40 x x 
41 vax X 
42 X 
43 X 
44 x xX 
45 x Xx 
4b x 
47 x 
1 2 3 4 


012345678901234567890123456789012345678901234567 


7 


WOON AU PWDY ~~ © 


Below is the chemical weapons attacks identifier matrix. 


1 2 5 4 
012345678901234567890123456789012345678901234567 
X X C C 


12 xX X X x 


~ 
lA 
< 
x 
4 


x 
x< 
Pad 


> 4 
x x 
=< 
=< 


=< 
=< 
> a4 
=< 
=< 


46 X 
47 X 


1 2 3 4 
012345678901234567890123456789012345678901234567 
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The matrix below shows which transition times were varied for the TBAD case given in Chapter V. 
(shown by the letter "T").The actual values for the TBAD case are given below the matrix. 


1 2 3 4 
012345678901234567890123456789012345678901234567 
X x X x 
X x X X 


OOBONAOUFWN — © 
< 
< 
< 
< 


—_ 
wi 
4 
= 

»< 


N 
— 
=< 
<< 
»~< 
< 


WW 
— 
< 
=< 
< 

x< 
»< 


46 X 
47 X 
1 2 3 4 


Ori-2 3 Meg 012s 4567890123456789012345678901234567 


= 
wi 
oO 


tGi2.13)= 
t(15, 15)= 
t(14,15)= 
EC1S 17)= 
t(16, 17)= 
t(17,29)= 
t(24,12)= 
t(25,13)= 
t(26,14)= 


165, t@h2, 14)= 
Soi, tC 15.20) = 
Poo, t(16,16)= 1.1, °t(16,17)= 0.935, t(14,26)= 0.55 
Boon Clo,2f)= 0.55 


0.55, t(12,15)= 0.825, t&12,24)= 0.275 
0 
1 
0 
.65, t(16,28)= 0.825, t(16,46)= 0.33 
0 
0 
0 
0 


“165 


pee (17 1h6)= 
55, t(24,25)= 
55, t(25,47)= 
55, t(26,27)= 


oo0oo}T-"00-"0=— 
— 
— 


Td 


t(27,15)= 0.55, t(27,47)= 0.77 
t(28,16)= 0.55, t(28,29)= 0.44 
t(29,17)= 0.55, t(29,47)= 0.77 


Wee 


The matrix below shows which transition times were varied for the TGOOD case given in Chapter 
V. (shown by the letter "T").The actual values for the TGOOD case are given below the matrix. 


1 2 3 4 
012345678901234567890123456789012345678901234567 
x x X x 


CON OA UE WN - oO 
> 
~< 
>o4 
> 4 


= 

W 

| 
=| 
+ 


—_ 
Wi 
=) 
— 

— 


46 X 
47 xX 
1 2 3 4 
012345678901234567890123456789012345678901234567 


Mle 15)= 450, tCi2,15)= 2.0, t(t2,24)= 0.75 
t(13,15)= 0.81, t(13.25)= 0.608, 1(13,47)= 0.7 
tose, 1o)= 420, t(14,17)= 2.0, t(14,26)= 1.5 
tase? )= i.0¢, t(l9,e7)= 1.22, 1(15,47)= 0.7 
t(16,17)= 4.0 

t(17,47)= 0.7 

t(24,12)= 1.0, t(24,25)= 6.0 

tice io)= 1.0, t(25,47)= 3.0 

t(2o-16)= 1.0, t(26,27)= 6.0 


19 


CG2(, to) = SO 6Ce7 (47) 0520 
t(28,16)= 1.0, t¢28,29)= 6.0 
t(29,17)= 1.0, t¢(29,47)= 3.0 
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APPENDIX G. APL FUNCTIONS USED IN THESIS 


This appendix contains the APL functions used for model 
analysis in this thesis. The APL used is the APL*PLUS System 
forethe PC from STSC, Inc. Some of the functions are 
documented and an introduction is included for each function 
in this appendix. Some of the functions are general in nature 
but the majority of them are "hardwired" to the specific model 
considered here and the author's particular hardware setup. 
Lines of code without bracketed line numbers are continuations 
of the lines above them. 

Thewt trot seunecLOnwAULDO-, 1s thegoverall function which 
solves the model and writes the output to the micro-computer's 
hard disk. AUTO uses the function FUND to solve for the 
fundamental matrix used in the solution technique from Taylor 


and Karlin. 
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ee 
E23 


C3] 


(4) 


CS) 


or! 

OTF 

ren 

Eo.) 

Pro 
C11] 
biZ 
Cis 
C14] 
fey 
C16] 
C17) 
alyssa) 
fare 
E203 
[233 
[22] 
E233 
C24] 
[25.1 
C26] 
C27] 


C28] 
P24 
1304 
ieee 
[327 
P57 
C34] 
354 


C36] 


OVR 'AUTO' 

T AUTO P3;A;0;R;CASE;SIZE;TABS; TSP; COMCLEAN ;s COMDIRTY ; DI 
SCLEAN ; DISDIRTY ; DECON; REFIT; LINE; LINE3 sy LINE5; LINE7:W;M 
USOJ ; LINED; OUTPUT 3;NUM;ALLTABS ; LT; OUTDET; TSPOUT;TSPMAT; 
C1:C2;C3;C4;:U;TABHEAD; TABNUM;TABOUT; DESCRIBE; FILE;HEAD 
ER;:N;ST;AT;TABLE;I1 


a THIS FUNCTION PROVIDES THE COMPLETE BREAKOUT OF WHER 
E A UNIT 

A SPENT ITS USEFUL LIFE. THE SOJOURN TIMES ARE COMBIN 
EDSINTO THE 

A FOLLOWING AREAS; COMBAT CLEAN, COMBAT DIRTY, DISENGA 
GED CLEAN, 

A DISENGAGED DIRTY, DECONTAMINATION AND REFIT. THE AC 
TUAL TIME AND 

A THE PERCENTAGE OF THE TOTAL UNIT LIFE ARE GIVEN. 


' INPUT THE NUMBER OF ABSORBING STATES ' 
Ac 


' INPUT THE CASE NAME FOR THIS RUN ' 
CASED 


A CREATE A FILE TO HOLD THE FUNCTION OUTPUT 


"INPUT THE NAME OF THE FILE FOR Your OUT Ur 

' THE NAMES CURRENTLY IN USE ARE SHOWN BELOW’ 
O<OFNAMES 

FILE¢D 

' “INPUT THE PILE TLE NUMBER VCR) OCh wr LL 

' THE TIE NUMBERS CURRENTLY IN USE ARE SHOWN BELOW’ 
O<OFNUMS 

NUM¢O 


FILE OFCREATE NUM 


A THE FIRST COMPONENT OF THE FILE IS THIS DESCRIBE VAR 
IABLE 
A WHICH DESCRIBES THE CONTENTS OF THESE TEE 


DESGCRIBEC 12 COG a 

DESCRIBE(CO;]¢€80T'THIS FILE CONTAINS THE ',CASE,’ RUN.’ 
DESCRIBE(1;J<80T'THE COMPONENTS IN THIS FILE ARE AS FO 
LLOWS: ' 

DESCRIBE(C(2;J]¢80T'COMP 1: THIS DESCRIPTION OF FILE CONT 
ENTS <7 

DESCRIBE(33;]<80T'COMP 2: THE ABSORPTION PROBABILITIES 

FOR THIS CASE.' 

DESCRIBE(4;]¢€80T'COMP 3: THE TIMES TO ABSORPTION ALL S 
TATES ' 

DESCRIBE[5;J]<80T'COMP 4: SUMMARIZED LIFE HISTORY OF UN 


82 


a7 


[38] 


£39] 


C40] 


C41] 


C42] 
C43] 
C44] 
C45] 
C46] 
C47] 
C48] 
C49] 
C50] 
tise Bee 
C32) 
3] 
[54] 
155 3] 
Esi6] 
ES? J 
ES6 J 
sees) 
C60] 
C61] 
eZ) 
C63 1 
C64] 
ces 1 
C66] 
C67] 
C68] 
C69] 
C70] 
Ey.) 
E72) 


[es J 
C74] 
E753) 
C76] 


E77) 
C78] 
yo 


A 


A 


A 


A 


A 


PrerOneeACH INITTA STATE.” 
DESCRIBE[6;]<80T’ THIS COMPONENT IS 
OF A 46xX15x75 CHARACTER’ 
DESGRIEEL7;I1<80T* MATRIX WHERE EACH 
ONDS TO AN INITIAL STATE. ' 
DESCRIEEL SO. 1CCOT COMPS 5750: 

UNIT FOR EACH INITIAL STATE.’ 
DESCRIBE[9;J]<¢<8s0T’ 

INITIAL STARTING STATE.’ 
DESCRIBEDIO; JesoTt’ BO GET THE REPORT FOR IN 
ITIAL STATE N REFER TO COMPONENT’ 

DESCRIBE cso N+ So 


IN THE FORM 
PAGE CORRESP 
DETAILED LIFE HISTORY OF 


EACH COMPONENT IS FOR AN 


OTCFF 
DESCRIBE ODFAPPEND NUM 

FIND THE FUNDAMENTAL MATRIX 
WeA FUND P 

PULL THE P MATRIX INTO ITS COMPONENT PARTS 
SIZE¢oP 
Q¢+(SIZE-A)TP 
R¢((SIZE[1]-A),0-A)TP 

ASSIGN THE PROBABILITIES OF ABSORPTION 
UeW+.xR 

FIND THE MEAN SOJOURN TIME FOR EACH STATE 
MUSOJ¢(0-A)1+/PxT 

THE TIME TO ABSORPTION FROM EACH INITIAL STATE 


ALLTABS¢ ,W+. X€ CoMUSOJ) ,1) eMUSOJ 
(C46 1 0(146)),UJOFAPPEND NUM 


TABHEAD¢ 6 70 0¢' ! 
TABHEAD[0;1¢70T'THE TOTAL TIME TO ABSORPTION FROM EACH 
OF THE INITIAL’ 


TABHEAD[1;]¢70T'STATES FOR THE ',CASE,’ CASE ' 
TABHEAD[2;]¢70T'IS SHOWN BELOW’ 

TABHEAD[33;]€700’ ' 

TABHEAD[4;]¢70T’' STATE TABS STATE 


TABS ' 
TABHEAD[5;1]€700’ ' 
ST¢R 2 23 9146 
AT¢& 2 23 eALLTABS 
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C80] 


rome) 
be2a 
Legg 
(84) 
C85] 
(86) 
C87) 
C88] 
C89] 
Bier 
(g2) 
[92] 
Eoey) 
C94] 
E954 
C96] 
Por) 
C98] 
Le) ) 
CE ileher) 
Galen 
[2024 
C103] 
C104) 
Dele) eet 
C106) 
aor 4 
eer e 
COI] 
C110] 
Cala 
Ca2 
Pas 3 


C114) 


Cio 
Creo 
C7 J 


bane?) 
CILILe) 3 
C203 
Ci24 
C122) 


Cue: 
C124) 
Paez) 


TABNUM¢'1I6,F15.4,110,F15.4’ OFMT(ST(;0);AT(;0);STL;:1); 
ATL 31) 
TABOUT¢TABHEAD, [0] 23 70 TTABNUM 


TABOUT OFAPPEND NUM 


A CREATE AN EMPTY BASKET TO HOLD THE GROUPED OUTPUT 
OUTPUT€ 46 15 70 oe’ ' 
INIT€0O 


TSPOUT+«( pMUSOJ) 0 
a START A LOOP 


LOOP1: 

A FIND THE TIME SPENT IN EACH STATE 
TSPeMUSOJXWLINIT; J 
A FIND THE TIMES FOR EACH CATEGORY 
COMCLEAN¢+/((601),4000)/TSP 
COMDIRTYe¢+/ (C1200) , (1291)5 220007 Tse 


DISCLEAN¢+/((690),(601),3400)/TSP 
DISDIRTY¢+/((2490),(12e1),1000)/TSP 


DECON¢+/( (3690), (601),490)/TSP 
REFIT¢+/((4200),491)/TSP 

A OUTPUT SECTION 

OUTPUTLINIT;0;]¢’ FOR THE ',CASE,’ CASE ',((7 
0-23)-pCASE)o’ ' 

OUTPUTLINIT;1;]¢’ AND INITIAL STATE OF ’,STA 


TNAMECINIT; J 
OUTPUTLINIT;2; 1¢€7 Oe) 


LINE3¢’ THE TOTAL UNIT LIFE UNTIL COMBAT INEF 
FECTIVE WAS ',(7 3 SALLTABSLINIT)),'’ HRS’ 
OUTPUTLINIT;3;]¢LINE3,(70-pLINE3)0’ ' 
OUTPUTLINIT;4;]¢700" ' 

LINE5¢’ OF ITS USEFUL LIFE THE UNIT SPENT IT 
IN THE FOLLOWING AREAS ' 
OUTPUTLINIT;5;]J¢LINE5,(70-pLINE5)o0’ ' 


OUTPUTLTCENIT ; 6 e7Oce 
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C126] 
CATED: 


C128] 
f129] 
C130] 
Cauca 
bis2 ) 


ia3.3 1 
C134] 
PEs 5 J 


C136] 
C137] 
C138] 


E239 J 
C140] 
C141] 


C142] 
C143] 
C144] 


C145] 
C146] 
C147] 


C148] 
C149] 
CabeyreR 
Galea ie 
C152) 
Pao. J 
C154] 
Clkspoye 
P56) 
C157] 
C158] 


C159] 


C160] 
[262] 
Calley 
C163] 
C164] 
C65) 


C166] A 


MEE) Gi h HOURS elas led 


NiwOF LIFE 


OUTPUT[ INIT;7;]¢LINE7,(70-oLINE7)oe' ' 

OUTPUTLINIT;8;1¢700" ' 

LINE«’ COMBAT CLEAN ',(6 3 SCOMCLEAN),’ 
', 6 2 6(100xCOMCLEAN+ALLTABSLINIT]) 

OUTPUTLINIT;9;]©LINE,(70-pLINE)p’ ' 

LINE¢’ COMBAT DIRTY ',(6 3 SCOMDIRTY),’ 
', 6 2 &(100xCOMDIRTY+ALLTABSLINIT]) 

OUTPUTLINIT;10;1]¢LINE,(70-pLINE)oe’ ' 

LINE¢’ DISENGAGED CLEAN ',(6 3 SDISCLEAN),’ 
', 6 2 6(100xDISCLEAN+ALLTABS[INIT)) 

OUTPUTLINIT;11;1]¢LINE,(70-oLINE)oe’ ' 

LINE¢’ DISENGAGED DIRTY ',(6 3 SDISDIRTY),’ 
', 6 2 6(100xDISDIRTY+ALLTABSLINIT]) 

OUTPUTLINIT;12;]«¢LINE,(70-pLINE)o’ ' 

LINE¢’ DECONTAMINATION ',.(6 3 SDECON),’ 

', 6 2 8(100xDECON+ALLTABSLINIT]) 

OUTPUTLINIT;13;]¢LINE,(70-pLINE)e’ ' 


REE Setoms GREPIT) . 
', 6 2 60C100XREFIT+ALLTABSCINITI) 
OUDEULEENIT; 14° 1¢LINE;, (C70-pLINEVe’ ' 


LINE¢’ 


TSPOuULel SPOUT ,TSP 


»>(C INIT=45) pENDOFLOOP1 
INIT¢CINIT+1 
>LOOP1 


ENDOFLOOP1:OUTPUT OFAPPEND NUM 


BECAUSE THE DETAILED OUTPUT IS TOO LARGE TO BE COMP 


UTED AS A MATRIX 


THIS SECTION COMPUTES AND FEEDS EACH PAGE OF OUTDET 
TO THE FILE 
AS A SEPARATE COMPONENT 
TSPMAT€ ( ( epMUSOJ) , p¢MUSOJ ) ep ( pPMUSOJ) JTSPOUT 
INIT¢O 


THIS SECTION PROVIDES THE MOST COMPLETE BREAKOUT OF 
WHERE A 
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[1e7) 


C168] 


C169] 
C170) 
C171] 
C172] 
C173] 
C174] 
hao 


C176) 
C7) 


C178] 


C1793 
C180) 
C613 


C182] 
C183) 
C184] 
C185] 
C186] 


(i872) 
Lise 
Pr39) 
C190] 
eel 


C192] 
Cale: sh 3 
C194) 


bo 5 3 
C196] 
Cio) 
[Se 2) 
Gib 
C200] 
C201) 
P2025 
C203] 
C204] 
C205) 
C206) 
C207] 
C208] 


A UNIT SPENT ITS USEFUL LIFE. EACH STATE IN WHICH TH 

E UNIT STAYED 

A FOR AT LEAST ONE PERCENT OF ITS LIFE IS LISTED ALON 

G WITH ITS 

A PROBABILITY OF ABSORPTION IN EACH FINAL STATE. 
HEADERG 10 Joc. = 

LOOP2: 

A THIS SECTION MAKES THE HEADERS FOR EACH PLANE OF TH 


E VARIABLE OUTDET 


HEADER[0;]¢’ 

=9 CASE) Oa 
HEADER 1; j<¢’ 
TNAMELCINIT; ] 
HEADER(2;J¢750° 


FOR "RHE" , CASE, "CASE, (( 75223) 


AND INITIAL STATE OF ',STA 


LINED¢' THE TOTAL UNIT LIFE UNTIL COMBAT INEF 
FECTIVE WAS ',(7 3 S@ALLTABSCINITJ),'’ HRS’ 
HEADER( 3; JeLINED, C7S-pLINEDI eG 


HEADER(4;1]€750' ' 
LINED¢’ OF ITS USEFUL LIFE THE UNIT SPENT 


IN THE FOLLOWING STATES ' 
HEADER(5; J¢LINED,(75-pLINED)o’ * 


LT, 


HEADER(UG; 1¢75¢6) ae 


GENE De. STATES 
2S TIME EACH TOTAL TIME’ 
HEADER(7;J¢LINED,(75-opLINED)Je’ ' 


VIS 


LINED¢’ 

VISIT IN STATE ' 
HEADER[8;]¢LINED,(75-pLINED)e’ ' 
HEADER[9;1¢75p' ' 


LT¢ , ALLTABS 

I«( (TSPMATCINIT; ]+LTLINIT])20.01) 
Ne+/I 

CleI4(~2 0)!STATNAME 
C2€(10 3)8((N,1)eCI/,WLINIT; ])) 
C3€(10 3)3(N,1)e(I/,MUSOJ) 

C4€(10 3)5((N,1)eC(I/,TSPMATLINIT;])) 
TABLE¢(N,75)1C1,C2,C3,C4 
OUTDET*HEADER, (0] TABLE 


OUTDET OFAPPEND NUM 
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209 J 

[210] »CINIT=45) ep ENDOFLOOP2 
Weed) =INITCINIT+1 

[212] »LOOP2 


23) 

(214] ENDOFLOOP2: 'THE COMPLETE RUN FOR THE ',CASE,' CASE IS 
DONE ' | 

(215] ° FOR A COMPLETE BREAKOUT OF WHAT THE UNIT DID SEE' 

C216] @eerheerilLE CALLED °,FILE.,’..' 

C27 ) 2 Ube E NUMBER SPOR SRAES Fie £5 ".,CSNUM),'’.' 

C218] PeOr ches DPSCKRIPTION OF THE FILESCONTENTS TYPE * O 
PREAD ‘9..0¢ (NUM, 1))5.0%.’ 

t229 ) 

V 


The following function, FUND, solves for the fundamental 


matrix W. This function is called by the function AUTO. 


DVR 'FUND' 
V WeA FUND P;30;31;3SIZE 

eT A THIS FUNCTION FINDS THE FUNDAMENTAL MATRIX GIVEN THE 
ONE STEP 

2] S PROBABILITY MATRIX, P, AND THE NUMER OF ABSORBING ST 
ATES, A, FOR 

C3] A A SEMI-MARKOV PROCESS. 

C4] A 

{5] a FIND THE SIZE OF THE Q MATRIX AND NEEDED IDENTITY MA 
TRIX 

C6] 

C7] SIZE¢CoP)J-A 

C8] Q+SIZETP 

[9] fhosiZhod,. (lLtolZb po 


[10] 
[11] a SOLVE FOR THE FUNDAMENTAL MATRIX W. 
[12] 
[13] WeB(I-Q) 
V 


The function CHEMPROB (following pages) is used to vary 


the probability of chemical agent attacks from various states. 
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C1] 
[23 


[3a 


C4] 
bod 
C6] 
Ca 
C8] 


[oo 


Eas] 


Fa] 


ba23 
p34 
C14] 
[aS 


[re] 


Gaby 


C18] 


C19] 
2a 
Ee 
P2200 


P2334 
[24] 
C25] 
C26] 
E27) 
P26 J 
Ze) 
bs04 
343 
[a2 
ies cya 
C34] 
L325) 
[360 
[3] 
P38) 


Vv 


A 


A 


A 


A 


VR  CHEMPROB: 
CHEMPROB P; FACTOR; B;SCHEM; DELTA; COUNT 


THIS FUNCTION ALLOWS THE USER TO VARY THE PROBABILIT 
OF 

COMING UNDER CHEMICAL ATTACK FOR ALL STATES THAT LEA 
Te 

POSSIBLE CHEM ATTACKS. 


INPUT SECTION 


' INPUT THE FACTOR TO CHANGE THE PROBABILITY OF CHEMI 
CAL ATTACK’ 

' FACTORS GREATER THAN ONE WILL INCREASE AND FACTORS 

LESS THAN ' 

' ONE WILL DECREASE THE OVERALL CHANCE OF CHEMICAL AT 
TACKS. : 

' INPUT FACTORS CAN BE BETWEEN O TO 3 AND DO NOT HAVE 
TOVSBE, INTEGER. 


FACTOR¢U 


' T0 VARY ONLY THE CHANCE OF CHEMICAL ATTACKS WHILE E 
NGAGED ' 

; INPUT *x*Ewun, 
ATTACKS ’' 

' WHILE NOT ENGAGED INPUT ***Crex, 
EOF 


' CHEMICAL ATTACKS FOR ALL POSSIBLE STATES INPUT ***A 
kkk! 


TO VARY THE CHANCE OF CHEMICAL ONLY 


TO VARY THE CHANC 


Be 


GUTS OF THE PROGRAM. 
A PARTICULAR 
STATE SPACE. 


THIS SECTION IS HARDWIRED FOR 


COUNT+¢0 
PNEW¢P 

JUMP TO THE CHEM ONLY SECTION IF REQUESTED 
+(B='C’ )@CHEMRESET 


LOOP TO DO THE ENGAGED CHEM ATTACKS 


ATTKLOOP: SCHEM¢€PL[ COUNT; (12+COUNT) 1+PL COUNT; (18+COUNT) J 


PNEW(E COUNT; (12+COUNT) J]€©FACTORxPL[ COUNT; (12+COUNT) J] 
PNEW[ COUNT; (€18+COUNT) J©FACTORXxP[ COUNT; (18+COUNT) ] 
DELTA*PNEW[ COUNT; (12+COUNT) ]+PNEWE COUNT; (18+COUNT)1-SC 
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[39 
C40] 
C41] 


C42] 


C43] 
[44] 
C45] 
[46] 
C47] 


C48] 
C49] 
[50] 
(aye 
C52] 
[53-3 
C54] 
(55) 
P64] 
Gay al 


[58] 


59) 


[60] 


[61] 
C62] 
C63] 
C64] 
C65] 
[66] 
E67) 
C68] 


C69] 
Ey O21 
bya 
[723 
C73] 
C74] 


[753 
[762 


HEM 
+ ( COUNT2=4 ) pBLOTTOCASE 


PNEWC COUNT ; €2+COUNT) J€ PL COUNT; (2+COUNT) ]-P[COUNT; (2+CO 
UNT) ]XDELTA+(1-SCHEM) 

PNEW[ COUNT; (6+COUNT) J€PLCOUNT; (6+COUNT) ]-PLCOUNT; (6+CO 
UNT ) ]x DELTA+(1-SCHEM) 

COUNT«<COUNT+1 

®ATTKLOOP 


BLOTTOCASE: 


PNEWL[ COUNT ; 46]¢PLCOUNT ; 46]-PL[ COUNT; 46]xDELTA=+(1-SCHEM) 
PNEW[ COUNT; €6+COUNT) J€ PLE COUNT; (6+COUNT) ]-PI COUNT; (6+CO 
UNT) ]xDELTA=(1-SCHEM) 

COUNT€<COUNT+1 

>+( COUNT<5) pATTKLOOP 


JUMP TO CONTINUE ***ALL*xs 
>(B='A' J oALL 


CALCULATIONS IF REQUESTED 


SECONDCASE: SCHEM€P[ COUNT ; (6+COUNT) 1+PLCOUNT; €12+COUNT) ] 


PNEW[E COUNT ; (6+COUNT ) J¢FACTORXxP[ COUNT ; (6+COUNT ) ] 

PHEW CCOUNT; C12+COUNT ) IGFACTORXPI[ COUNT; (12+COUNT) 1 
DELTA*¢PNEW[E COUNT; (6+COUNT ) ]+PNEW[ COUNT ; (12+COUNT) 1-SCH 
EM 

PNEWL COUNT; (COUNT-6) J€PL COUNT; (COUNT-6) 1-PLCOUNT; (COUN 
T-6) ]XDELTA+(C1-SCHEM) 

PNEWEC COUNT; (€18+COUNT) J€PLCOUNT; (€18+COUNT) ]-PLCOUNT; (18 
+COUNT ) ] XxDELTA+ (C1-SCHEM) 

PNEW[ COUNT ; (€24+COUNT) J€PL COUNT; (24+COUNT) ]-PLCOUNT; (24 
+COUNT ) ]X DELTA (1-SCHEM) 


COUNT+COUNT+1 
+ ( COUNT<7 ) pSECONDCASE 


THERDCASE: SCHEMCPL COUNT ; (6+COUNT) 1+P[£ COUNT; (12+COUNT) J 


PNEWL[ COUNT; (6+COUNT) J©FACTORXPLE COUNT ; (6+COUNT) J 

PNEWL[ COUNT; (12+COUNT) J€FACTORXPLE COUNT; (12+COUNT) ] 
DELTACPNEWL COUNT; (6+COUNT) 1+PNEWECOUNT; (12+COUNT) J-SCH 
EM 

FEW ECOUNT: (COUNT—6) 1¢PL COUNT; (COUNT-6)1]1-PL COUNT; CECOUN 
T-6) ]XDELTA+(1-SCHEM) 

PNEWLE COUNT; (€18+COUNT) J©PL[ COUNT; (18+COUNT) ]-PLCOUNT; (18 
+COUNT) ]XDELTA=+(C1-SCHEM) 

PNEWL[ COUNT; (24+COUNT) J€ PL COUNT; (24+COUNT) ]1-PLCCOUNT; (24 
+COUNT) ]XDELTA+(C1-SCHEM) 

PNEWC COUNT; (34+COUNT) J€P{ COUNT; (34+COUNT) ]-PLCOUNT; (34 
+COUNT ) ]xXxDELTA+ (1-SCHEM) 


COUNT«COUNT+1 


> (COUNT<11) epTHIRDCASE 
2END 
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ee 
C78] 
Gwe 
C80] 
C81] 
C82] 
ire oe 
C84] 
essed 


C86] 


C87] 


fee) 


feo4 
Cen J 
oa 
ioe 
eee] 
[94] 
Els) 


too 


ro7 4 


C98] 


Cet 1 


C100] 
ea oa 
C102] 
CPOs] 
C104] 
PasG5: J 


C106] 
Lae 7 1 


C108] 
eeoor 
Goce 
GI 
C1 a 


fies J 


C114] 


A 


ENTER INTO THE CHEM ONLY LOOP 


CHEMRESET: COUNT€6 


CHEM: SCHEM¢P[ COUNT; (18+COUNT) 1+PLCOUNT; (24+COUNT) J 


PNEW( COUNT; (18+COUNT) ]©FACTORxPL[ COUNT; (18+COUNT) ] 
PNEWL[ COUNT; (24+COUNT) ]©FACTORxP[ COUNT; (24+COUNT) ] 
DELTA¢ PNEWL COUNT ; (18+COUNT) ]+PNEW[ COUNT; (24+COUNT)1-SC 
HEM 

PNEWL COUNT; (COUNT-6) ]© PL COUNT; (COUNT-6 ) ]-PLCOUNT; (COUN 
T-6) ]xDELTA+(1-SCHEM) 

PNEWL COUNT; (6+COUNT) © PL COUNT; (6+COUNT) ]-PL£ COUNT; (6+CO 
UNT ) ]XDELTA+(1-SCHEM) 

PNEWL COUNT; (12+COUNT) ]©P[ COUNT; (12+COUNT) ]-PLCOUNT; (12 
+COUNT ) ]XDELTA+(1-SCHEM) 

COUNT+COUNT#+1 

+(COUNTs7 ) oCHEM 


CHEMREFIT: SCHEM¢P[ COUNT; (€18+COUNT) ]+P{COUNT ; (24+COUNT) ] 


PNEWC COUNT; (€18+COUNT) J©FACTORXPC COUNT; €18+COUNT ) J 
PNEWC COUNT; (24+COUNT) J©FACTORXP[ COUNT ; (24+COUNT ) ] 
DELTA¢PNEW[ COUNT; (18+COUNT) ]+PNEWE COUNT ; (24+COUNT)1-SC 
HEM 

PNEWCCOUNT; CCOUNT-6) JE P£LCOUNT ; C(COUNT-6) ]-PLCOUNT; (COUN 
T-6) ]XDELTA+(1-SCHEM) 

PNEWLE COUNT; (6+COUNT) J¢PLCOUNT; (6+COUNT) J-PLCOUNT; (6+CO 
UNT) J]XDELTA+ (1-SCHEM) 

PNEW([ COUNT; (€12+COUNT) J€PLCOUNT; €12+COUNT) J]-PLCOUNT; (12 
+COUNT) ] xDELTA+(1-SCHEM) 

PNEW[ COUNT ; (34+COUNT) J€PCCOUNT; €34+COUNT) J]-PLCOUNT; (34 
+COUNT) ] XxDELTA+ (1-SCHEM) 


COUNT¢COUNT+1 
> (CCOUNT<11) eCHEMREFIT 
>END 


A 
lave 


FINAL SECTION OF CALCULATIONS FOR THE ***ALL*** CAS 


ALL: SCHEM¢P( COUNT; (6+COUNT) ]+PL£ COUNT ; (12+COUNT) 1+PLCOU 
NT; €18+COUNT) 1+PCLCOUNT; (24+COUNT) J 
PNEWC COUNT; (6+COUNT) J€©FACTORXPL[ COUNT ; (6+COUNT ) J 
PNEW( COUNT ; (12+COUNT ) J©FACTORxP[ COUNT; (12+COUNT) J 
PNEWL COUNT ; (18+COUNT ) J€ FACTORXPL[ COUNT; (18+COUNT) J 
PNEW( COUNT ; (24+COUNT) J€FACTORXPL[ COUNT; (24+COUNT) J 
DELTA¢PNEW[ COUNT ; (6+COUNT) ]+PNEWLE COUNT; €12+COUNT) J+PN 
EWCCOUNT; (18+COUNT) ]+PNEW( COUNT; (24+COUNT ) ]-SCHEM 
PNEW([ COUNT ; (COUNT-6) J€ PLCOUNT ; (COUNT-6) ]-PLCOUNT; (COU 
NT-6) ] XDELTA+(1-SCHEM) 
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aS ') 
C116] 
Galen ara 
C118] 


aS] 
C120] 
[eZ 1 J 
tr22 I 
f¥23 J 


C124] 


C125] 


C126] 
27 | 
C128] 
CIS 
tas 0 J 


C131] 
rs 2 ) 
Vv 


COUNT+COUNT+1 
>» (COUNTS7 ) pALL 


ALLSEC: SCHEM¢P[ COUNT; (6+COUNT) ]1+PLCOUNT; (12+COUNT) 1+PLI 
COUNT; €18+COUNT) 1+PLCOUNT; (24+COUNT) ] 
PNEW[ COUNT ; (6+COUNT) J©FACTORxPI COUNT ; (6+COUNT) J 
PNEW[ COUNT ; (12+COUNT) J©FACTORXPLE COUNT; (12+COUNT) J 
PNEWL[ COUNT ; (18+COUNT) J©FACTORXxPLE COUNT; (€18+COUNT) J 
PNEWLE COUNT; (24+COUNT) J€FACTORXPL[ COUNT; (24+COUNT) J 
DELTA¢PNEW[ COUNT; (6+COUNT) 1+PNEWLC COUNT; €12+COUNT) 1+PN 
EWC COUNT ; (18+COUNT) ]+PNEWL COUNT; (24+COUNT) ]-SCHEM 
PNEWL[ COUNT ; (COUNT-6) J©PLCOUNT ; (COUNT-6) ]-PL COUNT; (COU 
NT-6) ]XDELTA+(1-SCHEM) 
PNEW[ COUNT ; (€34+COUNT) J©PL COUNT; (€34+COUNT) ]-PLCOUNT; (3 
4+COUNT) ]xDELTA+(C1-SCHEM) 


COUNT©COUNT+1 
+ (COUNT<s11)eALLSEC 


END: '’ THE PROGRAM IS DONE. 
BILITIES ARE IN’ 
' THE MATRIX PNEW' 


YOUR NEW TRANSITION PROBA 


The function CONPROB (following pages) is used in the same 


manner as CHEMPROB but is used to vary the probabilities for 


conventional weapons 


attacks. Both these functions are 


Pmreeractive. 


on 


OVR 'CONPROB'’ 
V CONPROB P;FACTOR;B;SUM; DELTA; C; ROW; NEWROW sy NEWPROB; DPRO 


B;ELSE 
(8) 
C2] a THIS FUNCTION ALLOWS THE USER TO VARY THE PROBABILIT 
Yeon 
[3] Aa COMING UNDER CONVENTIONAL ATTACK FOR ALL STATES THAT 
LEAD TO 
[4] aA POSSIBLE CONVENTIONAL WEAPONS ATTACKS. 


[5] 

C6) aA 

C7] 

C8] ' INPUT THE FACTOR TO CHANGE THE PROBABILITY OF CONVE 
NTIONAL ATTACK’ 

[9] ' FACTORS GREATER THAN ONE WILL INCREASE AND FACTORS 


INPUT SECTION 


LESS THAN ' 

C10] : ONE WILL DECREASE THE OVERALL CHANCE OF CONVENTIONA 
NL ATTACKS. ‘ 

baal ’ INPUT FACTORS CAN BE BETWEEN O TO 3 AND DO NOT HAVE 
TO BE INTEGER. ’ 

C12) 

[13] FACTOR¢O 

[14] 

eS al ’ TO VARY THE CHANCE OF CONVENTIONAL ONLY ATTACKS INP 
UT kakaCkkek, , 

[16] a: 

[awa '  T0O VARY THE CHANCE OF SIMULTANEOUS CONVENTIONAL’ 

C18] ' AND CHEMICAL ATTACKS INPUT *x«*Saxx, ' 

[19] oo 

C20] : TO VARY THE CHANCE OF CONVENTIONAL’ 

i234) : ATTACKS BOTH WITH AND WITHOUT CHEMICALS INPUT ***A* 
zk 

[22] 

[23] Be 

[24] 

[25] a GUTS OF THE PROGRAM. 

C26] A THIS SECTION USES THE METHOD FROM FUNCTION 

C27] A DP BECAUSE THERE IS ONLY ONE PROBABILITY PER ROW THA 
T NEEDS CHANGING. 

[28] 

[29] cCcé¢0 

C30] PNEW¢P 

[3A 

C32] 2%C€B='C’ )oCONVENONLY 

C33] 2>CB='S’ JoSIMUL 

[34] 2%ALL 

[353 

[36] CONVENONLY: 

C37] A PULL OUT THE ROW IN QUESTION 

[38] 

[39] CLEAN: ROW¢P[ (C+6);] 


o2 


C40] 
[41] 
C42] 
C43] 
[44] 
C45] 
[46] 
C47] 
[48] 
[49] 
C50] 
C51] 
52) 
Cee 
[54] 
fe> ) 
fio) 
[S73 
[58] 
feo) 
C60] 
rons 
C62] 
S627) 
C64] 
[65] 
C66] 


C67] 
C68] 
C69] 
C70] 
C71] 
C72] 
bigs: J 
C74] 
Wem 


C76] 


ea 7 | 


C78] 
f7od 
C80] 
C81] 
C82] 


C83] 
[84] 
C85] 


DPROBeP[ (C+6) 3:C) 
NEWPROB¢DPROBXFACTOR 


DELTA*¢NEWPROB-DPROB 

Puoee 1-DPROB 
NEWROW*ROW-ROWXDELTA+ELSE 
NEWROW[ CJ<¢NEWPROB 

PNEW[ (C+6) ; JENEWROW 


Cec+1 
+(C<5) oCLEAN 


DIRTY: ROW¢P[ (C+18) 35] 
DPROBeP[ (C+18) 3;C+6] 
NEWPROB¢CDPROBXFACTOR 


DELTA*NEWPROB- DPROB 
EESEG] -DPROB 
NEWROW*ROW-ROWXDELTA+ELSE 
NEWROW[C+6]¢«NEWPROB 

PNEW[ (C+18) ; J©NEWROW 


Cec+1 
4+(C<17)eDIRTY 
END 


A JUMP TO THE SIMULTANEOUS CHEM AND CONVENTIONAL ATTAC 
Koo OUnS TED 


Seu G sUMNe rt (C+6)3C+121+PL(C+6)>5Cr 1389 
PerPWieG+o6))-C+121¢FACTORXP[ (C+6)5C+12] 
PNEWICC+6)>C+18 JEPACTORXPI[ (C+6) 5C+18] 
DEDTAGPNEWE(C+6)3C+121+PNEW(E (C+6);C+18)]-SUM 


PHubWwece to). elePi(C+6)5C)—PlCet6) ;CIXxDELTA=C1-SUM) 
PNEWE (CC+6) 5€+241)¢€PL(C+6) 5;C+24]-PI(C+6) 5C+24]xDELTA+(1- 
SUM) 

Biwi €6+6)>C+307€Pi (C+6)):C+301—-Pi (C+6) :C+30)]xXDELTA=(1- 
SUM ) 

PNEW[ (€C+6) 5;C+40)¢PE (C+6) 5C+40]-PL (C+6) 3 C+40] XDELTA+(1- 
SUM) 

Seer]. 

7Lc=5)pSiIMUL 

END 


a DO THE CASE OF ALL CONVENTIONAL ATTACKS WHETHER COMB 
INED WITH 
A CHEMICALS OR NOT. 


MiieoUMNer i Ce+o)>7C)+P[ (C16) ,Cr123+Pi(C+6)>C+18) 
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PNEWL (€C+6)5C) ©PACTORXP( (Cro 

PNEW([ (C+6)5;C+121©¢FACTORXPI CExu6) Cai 

PNEW( (C+6);5C+18 JERACTORXPI (C¥G) erien 

DELTACPNEW[ (C+6);3;CJ+PNEW( (C+6)53C+12)+PNEWE (C+6) 5C+18I)- 


SUM 


PNEWL €C+6) 5;C+24 J) ¢€PC(C4+6)5C+24 J-PRCGrG) C+ 2a Liao 


SUM) 


PNEW( (€C+6);3C+30J¢PI (C+6) 5 C+30)-PL(C+6) Cr30j bers er 


SUM ) 


PNEW[ (C+6) ;C+40]¢€PL (C+6) 5C+40)]-PL (C+6) $C+40J]XDELTA+(1- 


SUM) 
CE€CrrE 


9(Cs5)eALL 


ALLTWO: ROW¢P[ (C+18) 35] 


DPROBGPICE+16)7Cr+6) 
NEWPROB¢DPROBXFACTOR 


DELTACNEWPROB-DPROB 
ELSE¢1-DPROB 
NEWROW¢ROW-ROWXDELTA+ELSE 
NEWROW[ C+6 J©NEWPROB 

PNEW[E €C+18) ; JENEWROW 


Cecis 


+(C<17) eALLTWO 


END 


END: 2 


THE PROGRAM IS DONE. 


BILITIES ARE 2h 
' THE MATRIX PNEW’ 


YOUR NEW TRANSITION PROBA 


The function DP (following page) allows the user to change 


one probability and still have the row of the P matrix sum to 


one. 


by the user remains constant. 
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The ratio between individual probabilities not targeted 


DVR 'DP' 
VY DPOUT<DP P;ROW; DPROB; NEWPROB; DELTA; ELSE; INDEX ; NEWROW 


A THIS IS AN INTERACTVIE FUNCTION THAT ALLOWS YOU TO C 
HANGE 

A THE PROBABILITY OF ENTERING DECON AT A LEVEL AND AUT 
OMAT ICALLY 

A UPDATE THE OTHER PROBABILITIES SO THAT THE ROW SUM R 
EMAINS ONE. 


A 


YOU WANT 
A TO CHANGE AS LONG AS YOU ONLY WANT TO CHANGE ONE PER 
A ae OF THE OTHER PROBABILITIES. 

DEOULCP 
START: ' INPUT THE ROW AND COLUMN OF THE PROBABILITY YO 


A 


U WANT TO CHANGE’ 


ACTUALLY THIS FUNCTION WILL WORK FOR ANY PROBABILITY 


INDEX<¢D 


ROW¢PCINDEX(0]; ] 
DPROB¢P[ INDEX(0); INDEX(1]] 


f 


DELTA¢NEWPROB- DPROB 
ELSECI-DPROS 
NEWROW*CROW-ROWXDELTA+ELSE 


PULL OUT THE ROW IN QUESTION 


THE OLD PROBABILITY IS 


', $DPROB 


ENTER THE NEW PROBABILITY’ 
NEWPROB¢O 


NEWROW[ INDEX(C1JJ<«NEWPROB 
DPOUTCINDEX([0]; J©NEWROW 


DO YOU WANT TO CHANGE ANOTHER PROBABILITY Y OR N’ 
9(O='Y’ JoSTART 


oS 


The functions CHEMRAND and CONRAND were used during the 


sensitivity analysis. These functions "protect" the 


probabilities of interest from the randomizing matrix used to 


perturb the rest of the P matrix. 
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The function REPORT is used to call up the Summary Report 


the initial State 6. 
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